Mode-Coupling Model of Mott Gap Collapse in the Cuprates: 
Natural Phase Boundary for Quantum Critical Points 
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A simple antiferromagnetic approach to the Mott transi- 
tion was recently shown to provide a satisfactory explanation 
for the Mott gap collapse with doping observed in photoe- 
mission experiments on electron-doped cuprates. Here this 
approach is extended in a number of ways. RPA, mode cou- 
pling (via self-consistent renormalization) , and (to a limited 
extent) self-consistent Born approximation calculations are 
compared to assess the roles of hot-spot fluctuations and in- 
teraction with spin waves. When fluctuations are included, 
the calculation satisfles the Mermin- Wagner theorem (Neel 
transition at T = only - unless interlayer coupling effects 
are included), and the mean-field gap and transition temper- 
ature are replaced by pseudogap and onset temperature. The 
model is in excellent agreement with experiments on the dop- 
ing dependence of both photoemission dispersion and mag- 
netic properties. The magnetic phase terminates in a quan- 
tum critical point (QCP), with a natural phase boundary for 
this QCP arising from hot-spot physics. 

Since the resulting T=0 antiferromagnetic transition is con- 
trolled by a generalized Stoner factor, an ansatz is made of 
dividing the Stoner factor up into a maieriaZ-dependent part, 
the bare susceptibility and a corre/afion-dependent part, the 
Hubbard U, which depends only weakly on doping. From the 
material dependent part of the interaction, it is possible to 
explain the striking differences between electron- and hole- 
doping, despite an approximate symmetry in the doping of 
the QCP. The slower divergence of the magnetic correlation 
length in hole doped cuprates may be an indication of more 
Mott-like physics. 

Discussions of interlayer coupling, doping dependence of 
U , extension to a three-band model, and polaronic effects are 
included. 



I. INTRODUCTION 

Schrieffer, Wen, and Zhang^ originally proposed that 
the magnetic insulating phase in underdoped cuprates 
could be understood via a spin density wave (SDW) 
approach to the Mott transition, and successfully de- 
scribed the spin wave spectrum of the undoped parent 
compound, which is an antiferromagnetic (AFM) insula- 
tor. Kampf and Schrieffer^ showed that precursors of 
the Mott transition could give rise to a pseudogap in 
the quasiparticle spectrum, between incipient upper and 
lower Hubbard bands (U/LHBs). Attempts were quickly 
made to go beyond mean field theories by incorporating 
fluctuation effects, but a number of problems soon arose. 
While some calculations found evidence for pseudogaps'^, 
others did not^. 



The rapid disappearence of Neel order with hole dop- 
ing created more problems: many calculations, even in- 
cluding strong fluctuations, predicted magnetic order as 
T — > 0, coupled with diverging magnetic correlation 
length ^, whereas ^ is found to remain finite even in 
the presence of the pseudogap. This has been used as 
evidence that the band structure picture of the Mott 
transition breaks down, and must be replaced by a local 
picture: 'Mott physics' instead of 'Slater physics'. On 
the other hand, other calculations find evidence for in- 
stabilities - either to incommensurate magnetism^ or to 
phase separation^' '', and the saturation of ^ could be due 
to nanoscale phase separation physics. The situation is 
at a stalemate, with some models neglecting both phase 
separation and magnetic effects, and explaining the pseu- 
dogap in terms of purely superconducting precursor ef- 
fects, while others find a magnetic quantum critical point 
(QCP) in the deeply underdoped regime, and yet others 
find a QCP above optimal doping. 

Clearly, a simpler alternative is an important desider- 
atum, and one has recently been proposed. While phase 
separation is a significant complication for hole doping^ 
this instability appears to be greatly reduced or absent 
in electron- doped materials^'^, allowing a much simpler 
analysis. Moreover, for electron doping, the band picture 
involving short-range commensurate AFM order seems 
justified, in that magnetic correlations remain commen- 
surate, while the correlation length diverges for all dop- 
ings up to the QCP. The desirability of a reference sys- 
tem free of phase separation complications coheres with 
Laughlin and Pines' observation^": "This problem [of 
identifying the correct quantum protectorate] is exacer- 
bated when the principles of self-organization ... com- 
pete. ... [Hjigher organizing principles are best identified 
in the limiting case in which the competition is turned off, 
and the key breakthroughs are almost always associated 
with the serendipitous discovery of such limits." 

While many models attempt to describe the properties 
of the cuprates over a limited doping range, it has proven 
difficult to systematically reproduce the changes over 
an extended doping range. Remarkably, simple mean 
field calculations^ were able to reproduce the full dop- 
ing dependence of ARPES spectra in the electron-doped 
cuprates^ in terms of a Mott gap collapse (QCP) near 
optimal doping. Here, these results are expanded upon 
in a number of ways. First, a number of models are ap- 
plied to the electron-doped system, to see the effects of 
various correlations. A key issue is finite temperature ef- 
fects: the RPA predicts a Neel temperature T/v ^ U - 
much larger than found experimentally. Proper inclusion 
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of thermal fluctuations, introduced via a self-consistent 
renormalization (SCR) model^^'^^, drives the Ncel tem- 
perature T^r to zero (Mermin- Wagner theorem), replac- 
ing the RPA gap Amf by a SCR pseudogap A* and the 

mean-field Neel temperature with a crossover tem- 
perature T*, with Amf ^ A*, ~ T*. Inclusion of 
spin wave scattering, via the self-consistent Born approxi- 
mation (SCBA) produces a large incoherent background, 
but the coherent part of the spectrum is recognizably 
the same as the spectrum found in the RPA and SCR 
approaches, with only moderate band renormalizations. 
The present conclusions (including Ref. 12) are consistent 
with more recent findings^^'^**. 

In summary, there is a QCP near optimal doping 
in the electron-doped cuprates, associated with Mott 
(pseudo)gap collapse. The transition is characterized 
by three concurrent factors: termination of a zero- 
temperature AFM transition (which can be associated 
with a finite-T Neel transition due to weak interlayer 
coupling); collapse of a pseudogap centered on (7r,7r); 
and crossover of the Fermi surface from small pockets 
to large barrel. Good agreement with experiment re- 
quires a weak Kanamori-style^^ renormalization of the 
Hubbard U with doping. The same model can describe 
both ARPES and magnetization results. A similar QCP 
is predicted at a comparable hole doping indeed a natu- 
ral phase boundary for magnetism exists, associated with 
hot spot physics. There is however, a striking difference 
in the hole doping case: saturation of a spin sum rule 
leads to much smaller correlation lengths and absence 
of finite-T Neel order. Stripe physics appears to play a 
lesser role - turning on at lower temperatures - possibly 
as a form of interaction of the doped polarons. 

This paper is organized as follows. Section II describes 
the SCR formalism. Related Appendices discuss the ex- 
tension to a three band model (Appendix A), the dop- 
ing dependence of U (Appendix B), and a more ac;curatc 
solution of the self-consistency equation (Appendix C). 
Since the transition occurs when a Stoner factor equals 
unity, it is controlled by the real part of the bare sus- 
ceptibility. Hence Section III reviews the properties of 
Rex, showing that plateaus in x as a function of doping, 
q, or Lu are all controlled by the physics of hot spots. In 
turn, these plateaus provide natural phase boundaries for 
QCPs. The resulting susceptibility has a form similar 
to that postulated for a nearly antiferromagnetic Fermi 
liquid (NAFL), but a calculation of the NAFL param- 
eters (Appendix D) finds that there are extra (cutoff) 
parameters, which cannot be neglected. In Section IV, 
this renormalized susceptibility is incorporated into the 
lowest-order correction to the electronic self energy, al- 
lowing a calculation of the spectral function associated 
with the pseudogap (Tjv = 0). Excellent agreement 
is found with the ARPES spectra of Nd2-i;Cej;Cu04±5 
(NCCO). A remaining problem lies in magnetic polaron 
effects which are expected at very low dopings: these 
bear some resemblance to nanoscale phase separation, 



and can lead to anomalous localization effects. Section 
V offers a brief introduction to these effects, by analyzing 
the self-consistent Born approximation at half filling. It is 
found that only minor quantitative changes to the earlier 
results are expected. An extension of the results to the 
hole-doped regime is considered in Section VI. The model 
also provides a good description of magnetic properties, 
as discussed in Section VII. Section VIII shows that inclu- 
sion of interlayer hopping leads to a finite (Appendix 
E). Results are discussed in Section IX, and Conclusions 
in Section X. Some of these results have been reported 
previously in the discussion of the mean-field results® and 
in a conference procedings^^. 

II. MODE-COUPLING CALCULATION 
A. Model Dispersion and Doping Dependence of U 

In the present paper the mean field results are ex- 
tended by incorporating fluctuations via mode-coupling 
theory following Moriya's self-consistent renormal- 
ization (SCR)^^'^'''^* procedure. Mode coupling theo- 
ries have been applied to charge density wave (CDW) 
systems^®'^", and have led to a successful theory of weak 
itinerant magnetic systoms^^'^^. They have also been 
used to study glass transitions^^, and recently extended 
to glasses in cuprates^^. The mode coupling analysis 
is particularly convenient, being the simplest model for 
which the Mermin- Wagner theorem is satisfied. The re- 
sulting pseudogaps compare well with recent photoemis- 
sion experiments in electron-doped cuprates. While the 
SCR technique can be generalized to deal with compet- 
ing phases^"^, only the antiferromagnetic fluctuations will 
be treated here. 

The cuprates are treated in a one-band model. By 
comparison with a 3-band model (Appendix A) , this can 
be shown to be an excellent approximation for the mag- 
netic properties. The bare electronic dispersion is 

ek = -2t{cx + Cy) - At'c^Cy, (1) 

with Cj = cos kiU. The dispersions for undoped 
Sr2Cu02Cl2 (SCOC) and electron-doped NCCO can be 
fit by assuming t = 0.326ey, t' /t = -0.276, with U 
taken as an effective doping dependent parameter®, with 
U = 6t at half filling. Similar parameters are found^^ 
to describe the spin wave spectrum^^ in La2Cu04: t = 
OMeV, t'/t = -0.25, and U/t = 6.2. The former values 
will be used here. 

Many textbooks on strong correlation physics^®'^^ note 
that the Hubbard U should be doping dependent, based 
on the original results of Kanamori^^, but there are no 
satisfactory results for the doping dependence in the 
cuprates. A simple model calculation, which gives semi- 
quantitative agreement with experiment in NCCO^'®, is 
described in Appendix B. 
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B. Self-Consistent Equation 

The SCR scheme is introduced to incorporate strong 

fluctuations near the antiferromagnctic wave vector Q. 
The (path integral) formalism is standard-^^ and only the 
main results are given here. The quartic Hubbard con- 
tribution to the Hamiltonian is decoupled by a Hubbard- 
Stratonovich transformation introducing spin wave fields 
(f). The Fermion fields are then integrated out, leaving 
an approximate quartic effective action, which describes 
fluctuations about the mean field solution due to mode 
coupling. In the SCR model, the dynamical susceptibility 
is found self-consistently as 



1 - Uxo{q,i(^n) + A' 
with the bare susceptibility 



e^^^ + uj + i5 



(2) 



(3) 



where 5 is a positive infinitesimal, and the RPA suscep- 
tibility given by Eq. 2 with A = 0. 

The leading divergence corresponds to AFM at g = Q, 
so the denominator of Eq. 2 - the (inverse) Stoner factor 
- is expanded in terms of the small parameters u) and 
q' = q — Q (analytically continuing iujn lo + ie): 



5q{uj) = 1 - Uxo{q-, w) + A = 5 + Aq'"^ - Blo"^ - iCu), 

(4) 

where 

d = l-Uxo{Q,0) + \ (5) 
and So = S — X. The self-consistent equation for 6 is 

^ = ^0 + ^ Do{q,iUn), (6) 

where w is a measure of the quartic mode-mode coupling 
(Appendix D4) and 



DQ\q,iuJn) = S + Aq'^ + C\u;n\. 



(7) 



The sum over Matsubara frequencies can be carried out 
using 



de 



ImX{e + i6) 



lUJn - e 



/■°° e e 
■ / d—coth—ImX{e + i6). 
Jo 7^ 2j 



(8) 



Then 



J2 Do{q,iu;n) 



(2^)2 



0V i 

q,lOJn 

w 2T (S + Aq'^r 



[Cef 



(9) 



Note the sharp energy cutoff in Eq. 9. This comes about 
because the linear-in-o; dissipation is a result of Landau 
damping of the spin waves by electrons near the hot 
spots, and therefore the dissipation cuts off when the 
spin wave spectrum gets out of the electron-hole contin- 
uum. The cutoff parameter is defined in Appendix 
D2, above Eq. DIO. Numerical calculations (Fig. 40) 
show that the cutoff can be quite sharp, particularly near 
the VHS. 



C. Approximate Solutions 

Equations 6, 9 can easily be solved in the limit T = 0. 
In this case, there is a transition at 



^0 = -12u 



dq'^a^ r-'^ de 





47r 



TT (Ag'2)2 



Ro = l-v, (10) 



= -ln[l + + 



tan ^(ttq) 



(11) 



with Uq = Aq^/a^. Since the right-hand side is finite 
and negative, fiuctuations reduce but in general do not 
ehminate the order at T = 0. At the RPA level (A = 0), 
the AFM instability is controlled by the Stoner crite- 
rion, ^0 — * 0. The quantum corrected Stoner criterion is 
Uxo ~ where representative values of rj are listed in 
Table I. 

However, for finite T, there are corrections ln{6), so 
6 cannot be set to zero, and there is no finite temperature 
transition (the Mermin- Wagner theorem is satisfied). To 
see this, it is adequate to approximate coth{x) as 1/x for 
X < 1 and 1 for a; > 1. In this case, Eq. 6 can be solved 
exactly. Appendix C. However, this exact solution is not 
very illuminating, and a simpler approximate solution 
will be given here. Since only the term proportional to T 
is singular, T and 6 can be set to zero in the remaining 
term. Defining 



So = So+T]-'^, 

Eq. 6 becomes 



(12) 



-—- / -^ian-i( ) 

2CT 



wA 



(13) 
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where the second line uses Eq. C6, below. Hence, there 
is no finite temperature phase transition, and 6 only ap- 
proaches zero asymptotically as T ^ 0: approximately. 



^ _ 2([7J'g-'r^l*o|/3uTo^ 



D. Susceptibility 



(14) 



Given the (inverse) Stoner factor 6q, Eq. 4, the 
renormalized susceptibility can be written in nearly- 
antiferromagnetic Fermi liquid (NAFL)^* form. 



x{q, ^) 



XQ 



with coefficients 
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(15) 



(16) 



(17) 



(18) 



(19) 



The similarity of Eq. 15 to the corresponding result for 
CDW's^" should be noted - the SCR is a form of mode 
coupling theory. 

In the renormalized classical regime, the vanishing of 
5 as T ^ is controlled by a correlation length, Eq. 17, 
which can be written as^^ 



27TpJkBT 



(20) 



Numerically solving Eq. 13 (or Eq. C7) for 6, then the 
spin stiffness Ps is exactly given by 



An 



(21) 



with ^0 = V Using Eq. 14, an approximate ps is: 



^1^0 1 



12?ia 



2 ■ 



(22) 



Ps is plotted in Fig. lib, with = 0.384ey, chosen 
to give a ps in agreement with experiment for x = 0, 
T — (Section VII). The T-dependcncc of the prefac- 
tor ^0 agrees with one-loop cr-model results'^" rather than 
the more accurate two-loop results^®'^^ . This difference 
is presumably a deficiency of the present model in not 
using fully self consistent parameters; it will be discussed 
further in Section VII. 



E. Parameter Evaluation 

The susceptibility Eq. 15 is well-known in NAFL^^'^^ 

and spin fcrmion'''^''''* theories and in rcnormaliza- 
tion group (RG) calculations of quantum phase 
transitions^^'^^. In these calculations, the parameters of 
Eq. 15 (equivalently, A, B, and C) are usually determined 
empirically from fits to experiments. However, the good 
agreement between experiment and mean field theory for 
electron doped cuprates encourages us to try to calculate 
these parameters from first principles, following Ref. 17, 
in terms of a renormalized Hubbard parameter Uef / and 
a mode coupling parameter u. This puts a special pre- 
mium on the bare susceptibility xo, which is assumed to 
control the main material, doping, and pressure depen- 
dence of the Mott transition, while the interaction pa- 
rameters U and u are relatively constant. In fact, it is 
found that U has a weak but important doping depen- 
dence, estimated in Appendix B, which is consistent with 
experiment. A single, doping-independent value of u is 
chosen to agree with t— J results at half filling. 

The motivation for this approach comes from experi- 
ence with another strongly correlated system: electron- 
hole droplets in photoexcited semiconductors. Here it 
was found^^'^* that the correlation effects were controlled 
by an isotropic density-dependent interaction potential, 
whereas the material, anisotropy, and uniaxial pressure 
dependence were controlled by the kinetic energy - i.e., 
by the bare band structure. A similar approach was ap- 
plied to CDW systems^^. 



III. HOT SPOT PLATEAUS AND GENERIC 
QCPS 

While the properties of xo are now reasonably well 
understood, they remarkably do not seem to have been 
used to derive the parameters of SCR or NAFL theory. 
Here this oversight is corrected. In particular, calcula- 
tion of the curvature parameter A is discussed below. A 
new cutoff parameter is introduced, which is essential 
in explaining the differences between the QCPs for hole 
and electron dopings. The corresponding frequency pa- 
rameter C, Eq D7, and its associated cutoff parameter 
a^j (below Eq DIO) are discussed in Appendix D2. 
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A. Plateaus in Doping Dependence 

1. Hot Spots 
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intersects the replica FS shifted by Q. They are located 
a,t Cx = —Cy = Cxo, with 



Cm = cos akxQ \ ttji 



(23) 



and equivalent points. The edges of the plateau are those 
points at which the overlap terminates (hot spots cease 
to exist). For the present band structure, hot spots exist 
only when the chemical potential fi is in the range At' < 
/i < 0, or for doping 0.25 > x > —0.19 (electron dopings 
are considered as negative). Since the two end points 
play an important role, it is convenient to label them, 
and they are here called 'hot' hot spot and 'cold' hot spot 
(or H-point and C-point) for the hole and electron-doped 
termination points, respectively. It will be demonstrated 
below that at each doping, the hot spots also lead to 
a susceptibility plateau in momentum space, around Q, 
collapsing to a logarithmic (square root) divergence at 
the H- (C-)point. The /f-point is the VHS, and hence 
also involves a conventional ETT. The physics is simpler 
near the C-point, where the topology hardly changes but 
the FS and Q-FS become decoupled (it is therefore a form 
of Kohn anomaly"'^). 



FIG. 1. (a) Susceptibility xo at Q as a function of doping 
for several temperatures. From highest to lowest curves near 
X = 0.1, the temperatures are T = 1, 100, 300, 600, 1000, 
2000, and 4000 K. Dotted line = 1/Ueff, dot-dashed line = 
1.5/Ueff. (b) Density of states Nf for the same temperatures, 
(c) Susceptibility xo at Q as a function of doping for several 
frequencies at T = IK: oj = 0.01, 0.1, 0.3, 0.6, 1.0 eV. (d) 
Pseudo-VHS (peak of xo) as a function of temperature Tv 
(circles) or scaled frequency = cj^/tt (squares); triangles 



In the self-consistent renormalization scheme, the T — 
AFM transition is controlled by a Stoner factor, 
C/i?e(xo) — where 77 > 1 includes a quantum correc- 
tion. Table I. Hence, the relevant quantity on which the 
study is based is the real part of the bare magnetic sus- 
ceptibility, Eq. 3. This susceptibility has been analyzed 
in a number of papers. Whereas usually only Im{x) is 
explored in detail (e.g., Refs. 40-42), Re{x) was studied 
in Ref. 43. The extended discussion which follows is in- 
tended to bring out salient features for the computation 
of the NAFL parameters. 

The doping dependence of xo{Qt'-^) is illustrated in 
Fig. la, where Q = (tTjTt). At low T, the susceptibility 
has a plateau shape, which is not present in the density 
of states, Np, Fig. lb. Beyond the plateau edges xo 
falls off sharply on both electron and hole doping sides 
of half filling. This sharp falloff explains the appearence 
of QCPs: the Stoner criterion is satisfied on the plateau, 
but fails when xo drops. 

The plateau shape is characteristic of hot spot physics. 
Hot spots are those points where the Fermi surface (FS) 



2. Mean Field Mott Transition 

For the parameter values expected in the cuprates, 
these susceptibility plateaus control the physics of the 
Mott gap collapse. As a function of doping, the mean 
field Mott gap is found to close at a doping just beyond 
the edge of the plateau, for both electron and hole dop- 
ing, Fig. 2. The solid and long dashed lines are the com- 
mensurate and incommensurate mean field Mott transi- 
tion temperatures T*{x) calculated using the estimated 
Ueff{x), dotted fine in Fig. 1. For electron doping, there 
is a double transition, first from commensurate to incom- 
mensurate antiferromagnetic order at the plateau edge, 
then to the loss of any magnetic order at a slightly higher 
doping (inset a). For hole doping, the dominant antifer- 
romagnetic order is incommensurate for all dopings, but 
the difference in Tjv becomes significant only near the H- 
point (inset b). When fluctuations are included (below), 
it is found that the Neel transition is shifted to zero tem- 
perature, while a pseudogap first appears near the mean 
field Tpf. Note that in the hole doped regime, there is 
good agreement between the mean field transition and 
the pseudogap (squares in Fig. lb = data of Krasnov^"', 
assuming 2A — 4.6T*). For the real cuprates, the ter- 
minations of the Mott gaps are preempted by supercon- 
ducting transitions, close to the critical regime. 
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FIG. 2. Mean field magnetic transition temperatures de- 
termined from Stoner criterion using Ueff of Fig. 1. Solid 
line: commensurate (at Q); long dashed line: incommensu- 
rate. Dot-dashed line = lOTjv , where Tjv is the onset of long 
range AFM order, from [ 45] and [ 46] (with filled circles). 
Insets = blowups near C- and H-points. Squares in inset b = 
pseudogap data of [ 44]. 

3. The Paeudo- VHS 

The susceptibility Fig. la has a remarkable doping de- 
pendence, with the large peak at the Van Hove singular- 
ity (VHS) shifting'*^ to half filling with increasing temper- 
ature T. The peak position of this 'pseudo-VHS' defines 
a temperature Tv{x)^ Fig. Id (circles). This behavior can 
readily be understood from the form of Xo(QjO), Eq. 3. 
The denominator — e^r^^j — ~4:t{cx + Cy), is indepen- 
dent oft', and hence has a stronger divergence than the 
density of states (dos). Indeed, this divergence matches 
the strong VHS found for t' = (perfect nesting), and 
like that VHS falls at half filling, x = 0. There is one cru- 
cial difference - at low temperatures, this divergence is 
cut off by the Fermi functions, which leave the integrand 
non zero in a wedge which intercepts the zone diagonal 
(where the denominator vanishes) only at isolated points: 
the hot spots. Hence, the residual divergence at low T 
is still dominated by the conventional VHS. However, at 
finite T, excitations along the zone diagonal become al- 
lowed, leading to a stronger divergence of Xo{Q,^) near 
X = 0. 

The strong temperature dependence of the pseudo- 
VHS is in strong contrast to the density of states, Np, 
Fig. lb, and also with the pairing correlations"'^. The de- 
nominator of the pairing susceptibility involves the sum 
of the energies, e-^ + e^r^g = —^t'cxCy, rather than their 
difference (as in Eq. 3), and hence always peaks at the 
ordinary VHS. 

The difference between nesting and pairing susceptibil- 
ities has a fundamental significance. By mixing electron 
and hole-like excitations, the superconducting gap is al- 
ways pinned to the Fermi level, and can open up a full 
gap at any doping. On the other hand, a nesting gap 



need not be centered on the Fermi surface, and is con- 
strained to obey Luttinger's theorem, conserving the net 
number of carriers in the resultant Fermi surface. Hence, 
the only way a nesting instability (such as antiferromag- 
netism) can open a full gap at the Fermi level is for the 
instability to migrate with increased coupling strength 
to integer filling of a superlattice zone (e.g., half filling of 
the normal state). 

Since the susceptibility has such a distinct temperature 
dependence from the density of states, one might ask how 
the frequency dependence compares. This is illustrated 
in Fig. Ic at low temperature (IK). While the frequency 
introduces additional sharp features and has an overall 
very distinct appearence from the T-dependence, never- 
theless the main peak also shifts from the VHS toward 
lower doping with increasing w - in fact, the shift is al- 
most the same when comparing fiLo and nksT, Fig. Id. 
The dashed line in Fig. Id is = hu;^ /nltB, with^" 

with T — 2t'/t and fi = ii/2t. The proportionality of 
frequency and temperature dependences holds only in the 
hole doped regime: temperature shifts the susceptibility 
peak only to half filling, a; = 0, while frequency will shift 
the peak beyond half filling (x < 0). 

The structure in the low temperature susceptibility. 
Fig. 1, with its largest peak at the H-point on the hole 
doped side, is in striking contrast to the calculated dop- 
ing dependence of the Neel transition. Fig. 2, which has 
a broad plateau on the electron-doped side, but falls off 
more quickly with hole doping, showing no sign of a peak 
near the VHS. This contrast can be accounted for by two 
effects. First, the shift of spectral weight with tempera- 
ture of the pseudo-VHS, noted in Fig. 1, would tend to 
produce a symmetric falloff of Tn with either electron 
or hole doping. But the dos peak at the VHS leads to 
better screening of Ueff for hole doping, thereby further 
depressing Tjs[. 

4- Neel Transition 

The mean field Neel transition is associated with short- 
range magnetic order, and hence should be compared 
to the experimental pseudogap transition T*, while the 
experimental Neel transition involves long-range mag- 
netic order. It is controlled by small parameters, such 
as anisotropy and interlayer coupling (Section VIII) and 
need have no connection to the mean field T/y. Never- 
theless, the mean field calculation provides an approxi- 
mate envelope of the resulting data, but overestimates 
the transition temperatures by a factor of 10, Fig. 2. 
The agreement is particularly good on the electron doped 
side (except for overestimating the doping of the QCP), 
while for hole doping the experimentally observed^^ 
(dot-dashed line) shows a stronger falloff, perhaps due 
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to phase separation. Since stripes can frustrate mag- 
netic order, the figure also includes the magnetic order- 
ing temperature of quasi-static stripe arrays, from Nd- 
substituted La2-2:Sr2:Cu04 (LSCO)^^, which is taken as 
a lower bound for the Nccl ordering transition in the ab- 
sence of stripes. A possible explanation for the rough 
proportionality of the mean field and long-range Neel 
transitions will be discussed in Section VIII. 



B. Plateaus in Momentum Space 

1. Plateaus 




FIG. 3. Susceptibility xo near Q for a variety of dopings at 
T = WOK. From highest to lowest solid curves near S = Q, 
the chemical potentials are fj, = -0.35, -0.30, -0.25, -0.20, -0.15, 
-0.10, -0.055, -0.02, and eV. For the dashed curves (top to 
bottom), /i = -0.352, -0.355, and -0.359eV. 

In analyzing either thermal fiuctuations or the quan- 
tum fiuctuations associated with QCPs, it is necessary to 
understand the susceptibility near the AFM vector Q. At 
each doping, hot spot physics leads to a plateau in mo- 
mentum space, centered on Q. Figure 3 shows how xo 
varies near Q at a low temperature (lOOK) for a series 
of different dopings. Results near T = are presented 
in Ref. 12. For all dopings there is a plateau in q. The 
width of the plateau at T = can be readily determined: 
in any direction, it is the minimum q needed to shift the 
replica FS so that the hot spots are eliminated. This 
can be found from the dispersion, Eq. 1, by substituting 
{Q + q)/2, or 

-2t{sx + Sy) - it'sxSy = fJ., (25) 

with Si — sm{qia/2). As shown in Fig. 4, this for- 
mula agrees with the (anisotropic) plateau width mea- 
sured from Fig. 3 (circles). The inset shows the shape 
of the plateau as a function of doping. The diamond 
shape of the plateau, Eq. 25, is related to the profile 
of the hole pockets formed by the overlap of the shifted 
and unshifted FSs. Specifically, the plateau is the re- 
gion of overlap of the two hole pockets, shifted to have 



a common center, as illustrated in Fig. 5. The remain- 
ing parts of the pockets also show up, as ridges''^ in the 
susceptibility, radiating from the corners of the diamond 
(similar to the peaks in the /i — 0.05ey data in Fig. 6, 
below). As noted by Benard, et al.^°, the susceptibil- 
ity in two-dimensions acts as a FS caliper. The plateau 
width leads to a natural limit on the magnetic correla- 
tion length, ~ I/Qc, in agreement with experimental 
data from YBaaCuaOr-a (YBCO)'^*'''^ (squares, trian- 
gles in Fig. 4), as noted previously^^'*^. Related data 
from LSCO^° are also shown. 
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FIG. 4. Plateau width q^, comparing Eq.25 (solid lines) 
and the measured widths (circles) from Fig. 3. Upper curve 
along [qc,0] direction, lower along [qc, 9c] /v^ direction. Sym- 
bols = experimental inverse correlation lengths from 
YBCO: large squares = Ref. 48, triangles = Ref. 49; LSCO: 
small squares = Ref. 50. Diamonds = TX/5000K . Dotted 
line: ^ = 100a. Inset = plateau boundary for a series of 
chemical potentials /i from (smallest) to -0.359eV (largest). 
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FIG. 5. Illustrating origin of plateaus (dotted line) from 
crossed hole pockets (short dashed lines). 
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2. Cusps 

For electron doping, the plateaus in q are particu- 
larly flat topped, Fig. 3. At low temperatures the edges 
sharpen up Fig. 6 and the fallofF in xa acquires a square- 
root singularity (Appendix Dl). The width of the plateau 
decreasing to zero as x ^ xc, and for electron-doping be- 
yond the C-point (/x > 0), the plateau ends and the sus- 
ceptibility displays split peaks away from Q, Fig. 6, with 
a dip in between. Thus xc is a QCP**^ where the mag- 
netic order changes from commensurate to incommensu- 
rate. (There is a corresponding QCP at the H-poinf*^.) 

However, the magnitude of xo also changes rapidly near 
/i = 0, so there should be an independent QCP from a 
magnetic to a non-magnetic phase near the same doping, 
as discussed in the previous subsection (note the line de- 
picting l/t/(/x = 0) in Fig. 6). 




(1.0,0.8) 



(1.0,1.0) (0.8,0.8) 



FIG. 6. Susceptibility xo near Q for several dopings near 
the C-point. Upper group a.t ^ = -O.OSeV, middle at ^ 
= (C-point), and bottom a.t fi — +0.05eV. Temperatures 
are T = 200K (dotted lines), lOOK (short dashed lines), 
lOK (long dashed lines), IK (solid lines). Horizontal line = 
U,ff{fi = Q). 

Technically, similar cusps also arise at the plateau 
edges for electron doping, > ^ > —0.22eV. The tops 
of the plateaus are not completely flat. Fig. 7a and the 
highest susceptibility is shifted away from Q (Appendix 
Dl). However, these effects are much weaker than those 
associated with ^ > {Ax/x ^ 0.5% - compare the 
vertical scales of Figs. 6, 7). Thus near the mean-field 
transition any structure on the plateaus is smeared out 
by thermal broadening. Even at T = 0, these features 
are likely to be negligible compared to dispersion in U 
which arises from renormalization effects^^. 



1.22 
(1 





FIG. 7. (a): Expanded view of susceptibility xo on the 
plateaus near Q for a variety of dopings at T = lOOK (solid 
curves^ or IK (dashed curves). From highest to lowest curves 
near Q, the chemical potentials are /i — -0.20, -0.15, and -0.05 
eV (for both solid and dashed curves). All curves except 
fi = —0.2061^ have been shifted vertically to fit within the 
expanded frame, (b): Similar plateaus for the hole doped 
materials (T = IK), with (from highest to lowest) fi — -0.359, 
-0.35, -0.3, -0.25, and -0.22 eV. 



3. Curvature (A) 

The plateau is a region of anomalously small lo- 
cal curvature A = A/U (Eq. 63) of the susceptibil- 
ity; Xo(Q + (?) = XQ ^ ^9^1 where A is an important 
NAFL parameter. Clearly, at T = lOOif the curvature 
A has gone negative near the iJ-point, Fig. 3. At even 
lower temperatures, it reverts to positive values. Fig. 7b. 
The temperature dependence of the normalized param- 
eter A' = {n I a)'^{A/t) is illustrated in Fig. 8 at several 
dopings. The temperature dependence is dominated by 
divergences at both H- and C-points. The divergence at 
the H-point Fig. 8a is the well-known logarithmic VHS. 
However, at finite temperatures spectral weight is shifted 
away from the VHS and A turns negative, only recovering 
a positive sign above T ~ 2{){){)K . The temperature at 
which A turns negative can be defined as Tincomm '■ A < 
for T > Tincomm- From Fig. Id, Ti„comm is comparable 
to but larger than Ty (for x < 0.06 A remains posi- 
tive). This in fact explains the origin of Tincomm- Fig- 
ure 8a demonstrates that A is negative at T ^ beyond 
the H-point (/i — —OAeV). Thus, increasing T above 
Tv produces the same susceptibility crossover. A similar 
crossover was discussed by Sachdev, et al.^^, except that 
they assumed that in the high temperature phase the 
AFM fluctuations remained centered on the commensu- 
rate Q, whereas here A is negative. At sufficiently high 
temperatures A again becomes positive for all dopings - 
i.e., the leading singularity of xo is always at Q. 

At the C-point, the collapse of the plateau width trans- 
lates into a divergence of the curvature at Q (A — s- oo). 
This divergence of the high-temperature susceptibility is 
cut off at low T, Fig. 8d, when the thermal smearing be- 
comes smaller than the plateau width. For smller T, A 
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is controlled by the curvature on the plateau. The tem- 
perature at which A has a peak, defined as T^, is plotted 
as diamonds in Fig, 4 (the peak is only found for a; < 0). 
Rather surprisingly, scales with the plateau width qc, 
even though the dynamic exponent is z = 2. Further, the 
maximum slope scales approximately as A„iax ^ T^~^'^ , 
which follows from the fact that A - T^^-^ at the C- 
point. 

At intermediate doping. Fig. 8b, c, A is generally a 
scaled-down version of the behavior near the two end 
points, with a crossover near fj, = —Q.25eV, where the 
T-dependence is weak. Also for intermediate tempera- 
tures, there can be fine structure on the plateau (e.g., 
solid lines in Fig. 7a) which can lead to wild swings in 
A{T). However, at these dopings they are not relevant, 
since the susceptibility peaks are away from Q, and this 
fine structure is not generally reported in Fig. 8. 




01 23401 234 



iog,.(T(K)) iog,„(T(K)) 

FIG. 8. Temperature dependence of A' for several dopings. 

C. Parameter Evaluation for Mode Coupling Theory 

The evaluation of the SCR parameters A and C was 
discussed above and in Appendix D. The collapse of the q 
and/or lu plateau widths near the H- and C-points leads 
to the introduction of additional parameters qc and a^. 
The narrow width of the g-plateau, particularly for elec- 
tron doping, leads to an additional complication not in- 
cluded in the conventional SCR analysis: the curvature 
of the bare susceptibility near Q — (tt, tt) (the S'-point 
of the BZ) is strongly temperature dependent, and for 
some dopings may even change sign. In principle, it is 
not difficult to incorporate an A{T) into the analysis near 
the mean- field Neel temperature (pseudogap onset). 
But for the present 2D system, long range Neel order 
only sets in at Tjv = 0, and for T << T^, a self consis- 
tent value of A should be found, by taking into account 
the effect of the pseudogap in modifying the electronic 
dispersion and hence x- For the present, this complica- 
tion is ignored, and in the following section A is taken 
as ^ = A{T^), where is the magnetic pseudogap on- 
set, the temperature where xo{Q)Ueff = 1, using the 
effective C/e// found earlier^ (Appendix B). This should 



be the most important A for controlling the pseudogap, 
and moreover at lower temperatures the band renormal- 
ization should strongly modify A{T). With this choice, 
the resulting A(/i) is plotted in Fig. 9a, along with the C 
parameter, evaluated at T = 0. For electron doping, this 
choice of A is always positive and varies smoothly with 
doping, diverging at the C-point. By contrast, for hole 
doping A is often negative, again illustrating the instabil- 
ity of the uniform AFM phase. Given A and C, Fig. 10 
shows the calculated values of Xq ^^^d ujsf, normalized 

to e- 




-10Lj I I I L_ 

-0.2 -0.1 0.0 „ 0.1 0.2 




FIG. 9. (a) Calculated values of A (circles for commen- 
surate Q, diamonds for incommensurate q ~ Q + q ') and C 
(squares). Solid line — Eq. D7. (b) Incommensurate wavevec- 
tor q' in two different directions: circles along (1,0), squares 
along (1, 1); dashed line: q' oc x. 

For hole doping, the incommensurability creates diffi- 
culties in defining the SCR model. The mean field transi- 
tion temperature at the incommensurate vector q is only 
marginally higher than that at Q, Fig. 2, suggesting that 
incommensurability should have only a small efect on the 
phase diagram. Thus, one might attempt to define a pos- 
itive A by maasuring the curvature from an incommensu- 
rate nesting vector. However, this A is highly anomalous, 
for a number of reasons. First, the incommensurate q' 
{q — Q + q ') forms roughly a square around Q, insert in 
Fig. 4, with the peak susceptibility generally along the 
(tt, 0) axis (circles in Fig. 9b). In this case, by symmetry 
there are four peaks in the susceptibility, at (tt ± q'^Tv) 
and at (tt, nzLq'). Moreover, the curvature measured from 
any incommensurate peak is highly anisotropic, since the 
susceptibility is nearly constant along the ridge of the 
square, with a shallow minimum at (7r,7r). Thus parallel 
to the ridge, A|| is nearly zero. Moreover, perpendicular 
to the ridge, A± takes on very different values on the 
sides of the ridge displaced toward or away from (7r,7r). 
The curvature is small, with significant deviations from 
quadratic (weaker curvature) moving toward (tt, tt), while 
moving away from (tt, tt), the curvature is larger, and de- 
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viating toward stronger curvature as the susceptibihty 
faUs off the edge of the plateau. For reference purposes, 
the average value of the quadratic part of A± is plotted 
as diamonds in Fig. 9a. (In this case, Lu^f has a peak 
near the VHS, Fig. 10.) This definition of A is almost 
certainly an overestimate. In the analysis of hole-doped 
cuprates in Sections VI and VII, the commensurate SCR 
model will be applied, with ^ as a free parameter. It 
will be found that agreement with measurement requires 
a somewhat smaller value for A than the estimated value 
of A±. This small A value, combined with the broad 
plateau, lead to a sum-rule saturation for xo smd a much 
slower divergence of £,{T) than found for electron doped 
cuprates. 
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FIG. 10. Calculated values of Xq/C^ (solid line) and cJs/C^ 
(short dashed line), assuming U = 6t. Long dashed line = 
doping x{fj,) (xlO); dot-dashed line = cosf^^ corrected for in- 
commensurate QT^Q (xl/5). 

In the following section, the present results are ap- 
plied to understanding the ARPES spectra of electron 
doped cuprates, concentrating on the four dopings an- 
alyzed by Armitage, et al.*. For convenience. Table I 
summarizes the parameters for these dopings. From the 
mean- field analyses^, the effective Hubbard parameters 
were found to be Uejf/t = 6 (a; = 0), 5 (a; = —0.04), 3 
(x = —0.10), and 2.5 {x = —0.15). [These numbers differ 
somewhat from those of Ref. 9, which included a second 
neighbor hopping, t" , to give the best fit of the Fermi 
surfaces.] The Stoner factor has a quantum correction, 
rj, Eq. 10, which tends to suppress the AFM transition; 
hence a smaller renormalization of U is required. This 
is reflected in Table 1: for x = -0.1, -0.15, there are two 
rows, the upper row using the mean-field U parameters, 
the lower with the quantum correction. Note that the 
L/'s are enhanced by essentially the quantum correction 
factor. These values will be used in the subsequent anal- 
ysis. 

The SCR analysis also involves a mode coupling pa- 
rameter u. An attempt to directly calculate u (Ap- 
pendix D4) failed, giving anomalously small values (Ta- 
ble I) due to the flatness of the susceptibility plateau 
{dx/dcu ~ 0), Fig. 39d. This problem has been noted 
previously, although there is debate about whether u 



diverges or vanishes^^'^^. Here ps is estimated from 
the measured correlation length for x — 0, using Eq. 20, 
as discussed in Section VII. Since ps oc u^^, Eq. 22, this 
gives = 0.384ey, which is assumed for all dopings. 
The calculated values of ps are illustrated in Fig. lib, 
based on Eqs. CI, 22. 



Table I: Electron Doped Cuprates 



X 


U/t 


A/a^ 


wi(eV) 






V 




(eV) 





6 


0.696 


0.345 


0.583 


0.635 


1.20 


1020 


760 


-0.04 


5 


1.16 


0.540 


0.455 


0.518 


1.17 


850 


3200 


-0.10 


3 


1.34 


1.32 


0.176 


0.342 


1.15 


500 


2700 




3.5 


1.56 


1.13 


0.206 




1.13 




2300 


-0.15 


2.5 


1.75 


2.16 


0.054 


0.172 


1.09 


56 


4000 


75 


2.9 


2.03 


1.86 


0.062 




1.05 


55 


3500 



It is convenient to compare the present results with 
parameters estimated for the SCR model^^ from experi- 
mental data for (optimally) hole-doped cuprates. The pa- 
rameters are defined as Tq = Aq'^/2nC, Ta = Aq-^ /2xo, 
yo = So{T = 0)/Aql, and yi ~ 12a^u/Tr^AC. The re- 
sults are listed in Table II, where the first line gives the 
hole-doped results estimated in Ref. 32. Moriya, et al.'^^ 
took q% = l/Ana^ (qsa = 0.282), while for Table II it is 
assumed that qs ^ Qc- A key difference is that Moriya, 
et al.'^^ assume the system is in the paramagnetic phase 
{yo > 0) at and above optimal (hole) doping, while in 
the present work yo < 0, and the system is paramagnetic 
due to the Mermin- Wagner theorem, with the Mott gap 
appearing as a pseudogap. The small magnitude of yo is 
suggestive of a system pinned close to a QCP. Finally, the 
parameter yi is estimated using the value = 0.384ey 
(above), and not the anomalous values of Table 1. 



Table II: SCR Parameters 



X 


To (K) 


Ta (K) 


yo 


yi 


- 0.2 


1600-4000 


3000-10000 


0.01-0.02 


3 


0.0 


180 


1150 


-5.27 


0.75 


-0.04 


310 


1300 


-3.31 


0.7 


-0.10 


380 


670 


-1.23 


1.5 


-0.15 


200 


220 


-0.31 


1.85 



IV. ARPES SPECTRA 

A. SCR Transition and Correlation Length 

Given the above parameters, the doping dependence 
of the MF and SCR transitions is compared in Fig. 11 
for the four electron dopings studied in Refs. 8, 9. The 
MF transition occurs when the bare Stoner factor So — 
1 — Xqo^ becomes negative, Fig. 11a. However, in SCR 
the renormalized Stoner factor S stays positive, so there 
is no T > phase transition (Mermin- Wagner theorem), 
although 6— So has a strong increase near the temperature 
where So changes sign. There is still a zero-T Neel transi- 
tion, controlled by the quantum corrected Stoner factor. 
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Sq = rj ~ XqqU ■ From Fig. 11c, it can be seen that at 
X — —0.15, the system is close to a QCP, 5q{T = 0) ^ 0. 
This QCP is controUed by the Stoner criterion of the 
zero-T antiferromagnet. While there is no long range 
order, there is still a Mott (pseudo)gap, controlled by 
short-range order, Fig. lid. A direct comparison of the 
transition temperatures is presented on a linear T scale 
in Fig. 12. The spin stiffness Ps (Fig. lib) is found to be 
nearly T-independent below the pseudogap onset. 
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FIG. 11. (a) 5— So (thin solid lines = — 5o); (b) ps calculated 
from Eqs. 22, 14; (c) -(5o; (d) A, Eq. 31. In all the plots, the 
solid curves correspond to x = 0.0, dotted lines: x — —0.04, 
short dashed lines: x = —0.10, long dashed lines: x = —0.15. 
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FIG. 12. Temperature dependence of gap A for (from high- 
est to lowest) a; = 0, -0.04, -0.10, and -0.15. Arrows show 
mean field transition temperature Tjv. 

While the value of U has been adjusted to fit the 
ARPES spectra, it is important to note that good agree- 
ment has also now been found with magnetic properties. 
This is discussed in Section VII. 



B. General Results 

Given the susceptibility, Eq. 15, the self energy can be 
calculated approximately as 



2 



E 



Ce 



with bare Green's function G'o(fc,ia;„) = l/(jw„ — ^j), 
= eg — and magnetic propagator Dq, Eq. 7; for 
the form of the integral, see the discussion near Eq. 9. 
In addition, xo = Xo(<3)0), q = Q + (f , n \s the Bose 
function, and 



1 ^ _ 3U 

(27) 



(Ref. 53). The last form is an approximation based on 
the empirical substitution xo 1/U in the pseudogap 
regime. [An improved approximation for E, (Go — > G in 
Eq. 26) is discussed in Section V.] After analytical con- 
tinuation, the imaginary part of the retarded self energy 
is 



ImY.^{k,u) 



V 



E 



de 



:ri(e) + /(4-+,-)]x 



Ge 



X'^(^ + ^-W (^ + ^g^)^ + (Ge) 



(28) 



The resulting self energy is plotted in Fig. 13 for T = 
lOOif . (The weak oscillations seen in some branches of 
S/ are an artifact due to an insufficient density of points 
in the numerical integration.) Note that /mS has the 
form of a broadened 5- function peaked at a; = ^k+Q' 



it were a (5-function, /mE = — ttA (5(a; — ^^r^g), then 



i?eS^(fc,tj) = - 

TT 



de- 



A" 



(29) 



so away from the (5-function 
1 



G(fc,cj) = 



a; - - i?el]«(fc, w) (^ - Cfe)(w " %+q) " 

(30) 



This is exactly the Green's function of the mean field 
calculation^ with the substitution A ^ A, where A 
can be evaluated by integrating 



A^ = 

TT 



dw/mS^(fc,w) 



(31) 



Fig. lid. This result is due to the Bose term n{e) in the 
square bracket of Eq. 28, the Fermi function / making 
no contribution. This leads to A being independent of k. 



(26) 
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FIG. 13. Imaginary part of the self energy, Eq. 28, assum- 
ing 1/C = 0.05t, 5 = 0.002, a„ = 1, T = WOK. The branches 
are labelled {kx, ky), in units of tt. 

Equations 30,31 constitute an important result, the 
connection between the Mott gap and short-range mag- 
netic order^^'^^. Recalling that A = J7 < Afj >, or 
A2 = [/2 < >2, where < M, >= (-1)* < 5"^ > is 
the staggered magnetization, then, in the spirit of an al- 
loy analogy, a short-range order parameter can be defined 
as 

-^0 <*,i> 

2 2 

= ^ E(^- + ''v)^+- ^ E 0) (32) 

which is equivalent to Eq. 31. (In the last equality in 
Eq. 32 the limit — > is an adiabatic approximation^^, 
while the approximation is made that x peaks near Q.) 
Thus, as long as there is short-range magnetic order (A. 
or ps non-zero), there will be a Mott (pseudo)gap. 

C. Application to the Cuprates 

Using the correct /mS^ from Eq. 28, and the calcu- 
lated parameter values from Table I, ARPES spectra are 
calculated for electron-doped cuprates, at the four dop- 
ings for which detailed data are available*. The resulting 
dispersions are shown in Fig. 14. There is a well defined 
pseudogap, with two peaks in the spectral function at a 
given k. It should be stressed that since there is no inter- 
layer coupling, long range antiferromagnetic order exists 
only at T = QK . The agreement with the mean field 
results^. Fig. 15, and experiment* is quite good, except 
that the SCR gap is smaller at half filling. This is due 
to lack of self-consistency: in calculating the self-energy, 
a susceptibility based on the bare Green's function was 
used, neglecting the opening of a gap near the Fermi 
level. In Section V it will be shown that when this is ac- 
counted for (via the self-consistent Born approximation) 
a larger gap is found. For completeness. Fig. 16 shows 



the mean field dispersion in the three-hand model, dis- 
cussed in Appendix A. The overall agreement in all cases 
is quite striking. 




FIG. 14. SCR Dispersion relations for electron doped ma- 
terials, calculated at T = lOOK: (a) a; = (U /t = 6), (b) 
X = -0.04 [U/t = 5), (c) X = -0.10 {U/t = 3.5), and (d) 
X — —0.15 {U/t — 2.9). Linewidth indicates relative inten- 
sity; for X — —0.15 all shadow features are extremely weak. 




FIG. 15. Mean field dispersion relations for electron doped 
materials, calculated at T = IK: (a) a; = {U/t = 6), (b) 
X = -0.04 {U/t = 5), (c) X = -0.10 {U/t = 3), and (d) 
X = —0.15 {U/t = 2.6). Linewidth indicates relative intensity. 




FIG. 16. Mean field dispersions in three-band model for 
electron doped materials, showing the two antibonding bands, 
assuming mg = 0.3 (a), 0.2 (b), 0.05 (c), and 0.01 (d). Other 
parameters are discussed in Appendix A. 

In an earlier calculation^^ a somewhat larger value of 
u was assumed, u^^ — 0.256eV. This leads to stronger 
quantum corrections: the parameter rj — 1 (Table I) was 
about twice as large and the gaps in Fig. 14 were smaller, 
particularly near half filling. 

Figure 17 shows typical calculated spectra for several 



12 



fc-points in the a-b plane. Broadened Hubbard bands are 
found, which gradually smear out at high temperatures 
as S increases decreases). 



(a) 



(b) 




FIG. 17. Spectral functions for (a) x = 0, (b) s = -0.04, 
(c) X = -0.10, and (d) x = -0.15, at T = lOOis:. Solid lines 
at (7r,0), and long dashed lines at (7r/2,7r/2). 

Figures 18- 20 illustrate the temperature dependence 
of Im{G) and Im{T,) for two dopings, x — and -0.15. 
The broadening of the peaks can be understood from 
Eq. 28: particle- hole excitations are present within a 
range zta^/C of C^+g- Away from this particle-hole con- 
tinuum the main peaks are sharp, while they broaden 
when they enter the continuum. 



IX. B). It is the divergence of this peak as T — > which 
signals the AFM transition. At low temperatures, the 
peak positions in Im{G) have a temperature dependence 
consistent with the collapse of the Mott gap - e.g., the 
LHB shifts to higher energies (toward midgap) at higher 
temperatures. Some experiments on hole doped cuprates 
find the opposite dependence^^, which can possibly be 
understood as a localization or phase separation effect. 
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FIG. 18. Temperature dependence of (a) spectral func- 
tion and (b) imaginary part of self energy, for x = 0.0 at 
(7r,0). Temperatures are 100, 500, 1000, 2000, 3000, 4000, 
and 5000K. 

Note that the Mott gap collapse is anisotropic: for 
the undoped case, the nodal gap collapses between 2- 
3000K, while a gap persists near (7r,0) above 5000K. 
/m(S) has striking oscillatory structure, particularly 
near (7r/2,7r/2), which produces a similar weak struc- 
ture in Im{G) at low T. [Similar, weaker oscillations are 
present near (tt, 0), which can be better seen in Fig. 4c 
of Ref. 12.] In addition, there is a very intense, strongly 
T-dependent peak in /m(E) exactly at Cfc+g- (Fig- 18b 
- also present but not shown in Fig. 19b - see Section 



FIG. 19. Temperature dependence of (a) spectral func- 
tion and (b) imaginary part of self energy, for x = 0.0 at 
(7r/2,7r/2). Temperatures are 100, 500, 1000, 2000, 3000, 
4000, and 5000K. 

In contrast, for x = —0.15, Fig. 19, the splittings are 
absent near {tt/2, 7r/2), and vanish near (vr, 0) by ~500K, 
and the lines actually sharpen on warming. If the effec- 
tive U is reduced to 2.5t, no splitting is found, but the 
peak position and broadening have an anomalous T de- 
pendence. Clearly, the system is very close to a QCP. 
Figure 21 shows in more detail how the spectrum evolves 
with U near this point. 
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FIG. 20. Temperature dependence of spectral function for 
X = -0.15 at (7r,0), for U /t = 2.9 (a) and 2.5 (c). Tem- 
peratures are 100 (solid line), 500 (long-dashed line), 1000 
(short-dashed line), and 2000K dot-dashed line), (b): imagi- 
nary part of self energy at T = lOOK, U/t = 2.9. 



FIG. 22. Fermi surface map for x = —0.10 (a) and —0.15 




0.20 



FIG. 21. 17-dependence of spectral functions for x = —0.15 
at T = lOOK near the T = QCP, for U/t = 2.5 (short dashed 
line), 2.7 (long dashed line), 2.9 (solid line), 3.0 (dot-dashed 
line), and 3.2 (dotted line). 

Finally, Fig. 22 displays Fermi surface maps for x = 
—0.10 and —0.15, showing the crossover from small to 
large Fermi surface. Hot spot effects are prominent at 
X = —0.15, pinning the Fermi surface to the zone diag- 
onal and broadening it at a pseudogap due to hot-spot 
scattering^^. These should be compared with the mean- 
field^ and experimental^ results. It should be noted that 
in the mean field calculation, it was necessary to include a 
t" parameter to reproduce the experimental hole pocket 
near the zone diagonal. Such a parameter would have 
shifted the Fermi surface across the zone diagonal, lead- 
ing to improved agreement with experiment here as well. 




0.0 0.2 0.4 0.6 0.8 1.0 
1.0 




(b). 



Thus, the SCR calculation agrees with the mean-field 
results®, if the mean- field gaps and transition tempera- 
tures are interpreted as the opening of a pseudogap at 
finite T, with the long-range AFM appearing only at 
T=0. Moreover, the overall dispersions, Fig. 14 are in 
quite good agreement with the mean field results® and 
experiments*. 



V. SELF-CONSISTENT BORN 
APPROXIMATION 

A limitation of the above calculations is that the self 
energy Eq. 26 is calculated using the bare susceptibil- 
ity, whereas the full susceptibility should be strongly 
modified by the opening of the Mott gap. This can be 
corrected for by including the full Green's function into 
the self energy calculation. This is conveniently done at 
T = 0, since there is a long-range ordered phase, and the 
sole difference between the SCR and mean-field calcula- 
tions is a weak renormalization of the Hubbard U. In 
this renormalized Green's function can be found 

by summing all the non-crossing diagrams; this is the 
self-consistent Born approximation (SCBA). By includ- 
ing the interaction of the quasiparticles with spin waves, 
it also incorporates the physics of magnetic polarons. 
Magnetic polarons closely resemble lattice polarons, lead- 
ing to both coherent and incoherent contributions to the 
spectral function, with considerable bandwidth renor- 
malization in the coherent spectrum. In the t — J model 
it is known that the SCBA gives a good description of 
exact diagonalization results on small lattices^®. 

Here, a simple calculation is presented to estimate the 
effect of the SCBA corrections on the dispersion of the 
insulating phase. Only the coherent band dispersion is 
included, and the SCBA is appHed to the RPA solution^°, 
which should be similar to the ordered SCR phase at T = 
0. The calculation of Chubukov and Morr^^ is extended 
to include both lower and upper Hubbard bands. The 
RPA dispersions of the upper (c) and lower (v) Hubbard 
bands can be written as follows. If the bare dispersion 
is ek ~ —2t{cos{kxa) + cos{kya)) — At'cos{kxa)cos{kya) — 

2t" {cos{2kxa) + cos{2kya)), then defining e\^^ — (e^ ± 
ek+Q)/2, Ei-^ = ^Je'^-^^ + A^, A = U<S,>, then 



El- = ei^^ ± Ei' 



(33) 



Here A is the AFM gap, A - U/2 at half filling. For 
large A, this can be expanded as 



Aqi cos kxttcoskya + 
-|-Ao2(cos 2fc2;a -|- cos2kya), 



(34) 



with 
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Aoi = J/2 ± t', 



(35) 



(44) 



A02 = J/2 ± 2t". 



(36) 



The same dispersion is found in the t — J model^^'^^^ 
suggesting that the SCBA will be an equally good ap- 
proximation here. 

The self-consistent equation (replacing Eq. 26) can be 
written 



6-1(^,0.) =c.-(£;r-M) 

^*c,i;(fc, q)G{k + q,OJ + OJq) 



(37) 



where ^c.« is a vertex correction for the upper (c) or 
lower (v) Hubbard band and Wg is the spin wave disper- 
sion. As will be seen below (Eq. 48), ^ oc t"^, so Eq. 37 
is independent of t, depending only on ratios t' /t, t" /t, 
and J/t. However, the final dispersion also scales with 
t, so any comparison with experiment requires all four 
parameters. For an arbitrary electronic dispersion, these 
quantities can be evaluated as follows. The transverse 
susceptibility in the RPA can be written as^'^^ 

x"'""(q,q,w) = 



x"^ (q,q + Q,'^) = 



Xq (q,t^) 



with 



r 1 1 



'-Ek + Ek+q - co Ek+ Ek+q + 



(40) 



and 



xh 



1 v^ A(4->+<)) 
2N ^ £;(-)£;(-) 

k ^k+q 
1 1 



^Ek+Ek+q-U) Ek+Ek+q+U)-' 

In the large U limit, Eq. 38 becomes 

1 1 , 



U! + Ulq Ul — Ulq 



with 




2 V aq+lq 



(41) 



(42) 



(43) 



7? 



cos{qxa) + cos{qya) 



(45) 



and 



J' J" 
aq = l + —{1- cos{qxa)cos{qya)) + — {1 - 72,), (46) 

with J = 4tVC/, J'/J = {t'/tf, J" I J = {t"/tf. As 
befits a Goldstone mode, = at g = (0, 0) and (tt, tt). 
In this case, '^cv = ^c,i;j with 

*c,. = - 4;i) ± ^,(4"^ + 4i). (47) 

with fiq = l/(2?7g), or 

^c,v = 16i [ ; ± (7fc - lk+q)\ ■ 



2 _ ^2 
q iQ 



(48) 



Given LUq, -E^'^, and "^c.v, Eq- 37 can be solved nu- 
merically to find both the coherent and incoherent parts 
of the ARPES spectral weight. However, the incoherent 

^llg^rt contributes to a weak background, and the exper- 
imental spectra arc generally compared to the coherent 
part. Hence, for present purposes what is needed is the 
dispersion of the coherent part of G. Following Chubukov 
and Morr^^ this can be simplified. The Green's function 
has the form 

(39) ^ 

G{k,u)) = = -, 

(49) 

with quasiparticle residue Z, band edge lOmax, damping 
7, dispersion Ek, with step fimction 0{x) = 1(0) for x 
> (<) 0. The quasiparticle residue can be found as 



1-Z 



d^q ^(fco,g) 

47r2 {Ojq + Ek„+q)^ ' 



(50) 



where ko is the band-edge momentum: ko = {-k 12,11/2) 
[(tt, 0)] for hole [electron] doping. (With the conventional 
signs t' < 0, t" > 0; in the special case t' = t" ~ 0, both 
energies are degenerate.) An equation for the dispersion 
can then be found by substituting Eq. 49 into Eq. 37, 
and setting ui = iVmax- 



2„c,v 



^c,v ^ f fq_r_Hk^q)_ _ ^(fco,?) 

J 47r2 ^UJq + Ek+q UJq + Ek„+q 



(51) 



]], (52) 



(The damping adds a small correction to the dispersion, 
which we ignore.) It is convenient to rewrite Eq. 34 as 
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^k '" ~ 4^01 (coi + cos{kxa)cos{kya)) 
+^02(co2 + cos{2kxa) + cos{2kya)), 



(53) 



with coi = (1), co2 = 2 (-2) for the fower (upper) Hub- 
bard band. It is found that Ek satisfies a similar equa- 
tion, with renormalized Aoi Ai. In this case, the self- 
consistent equation Eq. 51 can be reduced to a pair of 
equations at fixed fc-values. For example, at fc = (0, 0) 

i{l + coi)Ai + (2 + CQ2)A2 
= Z[4(l + coi)Aoi + {2 + co2)Ao2] - Z^e'^^^^^ (54) 

with a similar equation at k — (tt, 0) [or (7r/2, 7r/2)]. Fig- 
ure 23 illustrates the self-consistent values of Z, Ai, and 
A2 as a function of J for fixed t' , t" . 
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FIG. 23. Renormalized parameters for lower (a) and up- 
per (b) Hubbard bands: Z (solid lines), A\ (dot-dashed 
lines), and A-i (dotted lines), in comparison with ZAi)\ 
(long-dashed-short-dashed lines), and ^Ao2 (short-dashed 
lines), assuming parameters t — 0.326eV^, t' — — 0.375f, and 
t" = 0.15t. Horizontal lines = experimental range for Ai 
(long dashed lines) and A2 (solid lines), after Refs. [ 64, 
65] . Also shown are the individual renormalization factors Zi 
(long-dash-dotted line) and Z2 (long-dash-dot-dotted line). 
Vertical lines delimit parameter values consistent with exper- 
iment. 

Note that any attempt to extract the bare parameters 
from the measured dispersion is highly underdetermined. 
Thus, while the band dispersion Ek and spin wave dis- 
persion LUq depend explicitly on J, t' , and i", the ver- 
tex function depends on t, so there are four parameters 
to determine, but only two parameters Ai and A2 can 
be found from the ARPES dispersion. Moreover, from 
Fig. 23a, the value Ai is insensitive to J in the range of 
interest. In principle, the parameters can be determined 
from additional measurements, including the Mott gap 
A, the spin wave velocity (as q ^ 0, ujg ^ Cgq), 



= 2a 



J(^ + J' + 2J"), 



(55) 



or the maxima in the spin wave spectra, i^(7r/2 7r/2) = 
2( J J' 2J"), W(^,o) = 2(J 2J'). 

Given this indeterminancy, a simplified picture is as- 
sumed here to estimate parameter changes: the renor- 
malized value oft — 0.326ey is assumed fixed, to keep ^' 
and the experimental ratios Ai/t constant, and further, 
the ratio t' /t" = —2.5 was assumed constant. Then the 
pairs of solid and long-dashed horizontal lines in Fig. 23a 
give the experimental ranges^"''^^ for Ai and A2 respec- 
tively. A reasonable match can be found for a bare 
t' ~ — 0.375t. In this case, the value of A2 suggests a 
bare J in the range 0.33 - 0.41i, or 108 - 135mey. For 
the same parameter range, the individual parameters are 
renormalized by = A^/Aq,, with Zi ~ 0.059-0.020 [0.8- 
0.84], Z2 ~ 0.44-0.51 [-0.21 - -0.54] for the lower [upper] 
Hubbard band. The ratio j/t" must be renormalized by 
the SCBA, since A2 and ^02 cross zero at different values 
of J, causing Z2 to be negative for the upper Hubbard 
band (it diverges when A02 0). 
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FIG. 24. Parameter values consistent with ARPES data 
for J (a) and 2 A (b). Broad range determined by Ai; narrow 
range (in both a and b) by A2. Circle — SCR result, square — 
best SCBA approximation. Dashed line in a: J = 0.6t—1.5t" . 

The above calculation can be repeated for different val- 
ues of t", and the allowed parameter values for J and t" 
are shown in Fig. 24a. The 'best' SCBA value (square) 
differs from the SCR value (circle) by less than a factor of 
two. Since J is reduced by polaron coupling, U = it'^/J 
must increase, Fig. 24b. This can be seen directly from 
the self-consistent equation for G. The leading edge of 
the band is found from Re{G~^{kQ,ujmax)) = 0, or 



47r2 



(56) 



The gap 2A is equal to the splitting between the upper 
and lower Hubbard bands at (7r/2, 7r/2) - it is noi the sum 
of the Wmas's for these two bands, since the bottom of 
the upper Hubbard band lies at (tt, 0). Correcting for the 
renormalization of the dispersion at (7r/2,7r/2) reduces 
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the gap, but even so the renormahzed A (Fig. 25, short- 
dashed hne) is larger than the bare value (solid line). 
Figure 26 shows that the SCBA increases the Mott gap 
near half filling, which corrects a shortcoming of the SCR 
model, noted above. 




FIG. 25. Renormalized U parameter as a function of J: 
solid line = bare U — it^ /J in large gap limit; short-dashed 
line — renormalized U from Eq. 56; long-dashed line = bare 
A corrected for the small gap limit, Eq. 56; dotted line = 
renormalized U in the small gap limit from Eq. 56; dot-dashed 
line = Ao3. 
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FIG. 26. Comparison of mean-field (dashed lines) and 
SCBA (solid lines) dispersions for t' = -0.375t, t" = 0.15t, 
and two choices of J, J/t — 0.42 (a) or 0.33 (b). 

A. Extension to Small U 

The above results were valid for the large-[/ limit, 
where the gap parameter A = U < Szi >^ U/2 >> t, 
in which case J — 2t'^/ A. As A decreases, certain modi- 
fications are necessary. The most important is a modifi- 
cation of J. From the above analysis, the susceptibility, 
spin wave dispersion, and renormalized band parameters 
all depended on the bare electronic dispersion. Hence, 
the value of J should be chosen to best approximate 
the bare Api, Eqs. 35, 36. This can be accomplished by 
matching the exact dispersion to the approximate form 
at fc = (0,0), or 



Fig. 27. Note that J ^ 1 as A ^ 0, Fig. 25. It is inter- 
esting to note that when the renormalization correction, 
Eq. 56, is added in, the renormalized A (dotted line in 
Fig. 25) lies close to the perturbative result A — 2t^/J 
(solid line). 
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FIG. 27. Mean field band structure (solid lines) plus ap- 
proximations involving Eqs. 57, 58, with (short-dashed lines) 
or without (long-dashed lines) a finite Aoa, for A/t = 0.5 (a), 
1.0 (b), 2.0 (c), and 3.0 (d). 

For small A an additional correction is required, to ac- 
count for quartic corrections in t/A. This can be done, 
as above, by adding a term ^93(1 — C2xC2y) to the model 
bare dispersion, which allows a fit to the exact bare dis- 
persion at (7r/2,0), if 

This yields a very good approximation to the dispersion 
down to A = t/2, Fig. 27 (short dashed line). The pa- 
rameter is plotted in Fig. 25. 

B. Summary of SCBA Results 

(1) Thus at half filling polaronic effects renormalize the 
bandwidth by only a factor of ~2, with some change in 
lineshape. Polaronic effects reduce the values of J, t' , and 
t" , and hence increase the value of U. Thus the gap is 
enhanced at half filling, correcting a shortcoming of the 
SCR calculation (Figs. 14, 15). 
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While the present calculations are restricted to half fill- 
ing, some additional features can be extracted from the 
calculations of Kusunose and Rice (KR)^'^. (2) While 
considerable weight is transferred to an incoherent spec- 
trum, the coherent spectrum is quite similar to that 
found in RPA and SCR calculations, and it is this com- 
ponent which is mainly seen in the ARPES spectra. Pos- 
sible evidence for the incoherent states is a second peak 
seen in ARPES spectra of half-filled cuprates^^, about 
0.6eV below the main peak of the LHB near the nodal 
point. While KR find an incoherent peak at half fill- 
ing about 2.5t below the first peak, its intensity actually 
maximizes away from the nodal direction toward F, while 
the experimental peak is stronger in the opposite direc- 
tion, towards {tt,tt). 

(3) An important result of the RPA and SCR calcula- 
tions is that U must decrease with doping to reproduce 
the experimentally observed crossover to a large Fermi 
surface. The same result has been found by Senechal 
and Tremblay^''. The results of KR are consistent, in 
that KR kept U doping independent, and did not find 
this crossover. 

(4) Whereas in the mean-field and SCR calculations, 
electron doping shifts the Fermi level into the UHB with- 
out affecting the relative weights of the two subbands, 
in the SCBA the UHB states below the Fermi level are 
formed by spectral weight transfer from the LHB. This 
spectral weight transfer had been seen experimentally, 
and its absence was known to be a shortcoming of mean 
field theory, which is thus seen to be corrected in the 
SCBA. 

(4) In hghtly doped NCCO, Armitage, et al.* found 
an additional weak pseudogap - actually a leading edge 
gap at the Fermi level of the UHB - which was not repro- 
duced by the SCR calculation. Such a pseudogap is found 
by KR, and in an earlier calculation by Stanescu and 
Phillips^"^. KR interpreted this as evidence that the filled 
states were not actually part of the UHB but were in-gap 
states close to the bottom of the UHB. Similar in-gap 
states had been proposed for hole-doped cuprates^^, and 
have been considered as evidence for stripes. [In LSCO, 
where stripes are most clearly observed, the added states 
are close to mid-gap^^; in other hole-doped cuprates, the 
evidence is less clear, but if in-gap states exist, they must 
lie close to the top of the LHB.] 

The connection between polarons and stripes is a del- 
icate issue: for very light doping one would expect 
magnetic polarons to form for both hole and electron 
doping. These polarons are strongly dressed electrons, 
with many features of second-phase inclusions, and have 
been suggested to act as precursors for nanoscale phase 
separation^°'^^. In hole-doped cuprates, there is consid- 
erable evidence that these polarons tend to cluster and 
form stripes. In electron-doped cuprates there is con- 
siderably less evidence for stripes, and it may be that 
polarons do not form clusters. Hence, the differences 
between a polaronic phase and a stripe phase might be 
rather subtle. 



VI. EXTENSION TO HOLE DOPED CUPRATES 

Thus, for electron-doped cuprates, a three-fold coin- 
cidence of Mott gap collapse, Fermi surface crossover, 
and zero-T QCP is found. SCR theory predicts a similar 
triad for the hole doped cuprates, and the present section 
explores the extent to which this is found experimentally. 

A. Pseudogap 

In hole doped cuprates, ARPES finds two features 
which are commonly referred to as pseudogaps ~ a 'hump' 
feature found near (tt, 0) at higher binding energy than 
the main, superconducting 'peak', and the 'leading edge 
gap', a loss of spectral weight in the immediate vicin- 
ity of the Fermi level. This latter feature is not ex- 
plained by the present calculation; it may be the mag- 
netic feature discussed in Section V.B^''''^'^'^'', or it may 
be associated with the onset of strong superconducting 
fiuctuations^^'''^. 

On the other hand, the 'hump' feature can be 
consistently interpreted as the collapse of the Mott 
pseudogap^^. Bilayer splitting cannot explain SIN tun- 
neling measurements^^'^^''*^ which find two hump-like 
features, roughly symmetric about the Fermi level. Cor- 
relation with ARPES suggests that the tunneling peaks 
refiect structure near (tt, 0), and Figure 28 shows that 
semi-quantitative agreement with experiment can be at- 
tained in terms of weakly split Hubbard bands, for a 
screened U — 2.3t (see also Fig. 2b). For simplicity, the 
calculation is carried out at the mean-field level. Fig- 
ure 28c, d shows how the bottom of the UHB near (tt, 0) 
gradually merges into the VHS of the LHB. The inten- 
sities and positions of the two dos peaks reveal a clear 
asymmetry. As the Mott gap vanishes, the two peaks 
merge into the VHS of the bare band. (There may be 
complications due to nanoscale phase segregation, since 
STM studies suggest that the peak and hump features 
are spatially segregated^^.) 
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FIG. 28. (a) Calculated dos for a series of hole doped 
cuprates, Eissuming Ueff = 2.3t, with x = 0.176 (solid line), 
0.184 (long dashed line), 0.202 (short dashed line), 0.225 (dot- 
ted line), and 0.244 (dot-dashod lino), (b) Comparison of shift 
of lower dos peak (circles) from (a) with representative tun- 
neling data [ 44] (triangles). (c,d) Band dispersion near the 
pseudogap for x = 0.176 (c) and 0.244 (d). 

The above interpretation requires that for hole doping 
also the Mott gap must collapse slightly above optimal 
doping. This is consistent with recent experimental ob- 
servations of a QCP''^. Moreover, the model predicts 
that at the QCP, where the pseudogap just closes, the 
Fermi level is exactly at the VHS (H-point). This result 
had been found experimentally in some lightly overdoped 
cuprates''^' 

In a recent confirmation of the QCP just beyond opti- 
mal hole doping''*, it is suggested that this 'has to do with 
the restoration of the Fermi-liquid state in the overdoped 
regime characterized by a large Fermi surface' (emphasis 
added); a similar conclusion was made by Balakirev, et 
al.^^. Thus two elements of the QCP triad are present. 
The third is more elusive. 



B. T=0 QCP 

In electron-doped cuprates, a finite Ncel temperature 
persists all the way to the QCP; by contrast, for hole 
doping Tjv — > at a doping x 0.02 — 0.03, considerably 
below the proposed QCP. Here it is suggested that a T=0 
magnetic transition persists out to the QCP, but the cor- 
relation length grows so slowly that three-dimensional 
Neel order is superceded by the superconducting tran- 
sition. Details are presented in a related publication^^, 
and only briefly summarized here. 

The key insight is that the susceptibility must satisfy 
the fluctuation-dissipation theorem^^'^^ 



with a and b constants. From Eqs. 6, 9, it can be shown 
that 



<-^>=-/^h/(0(^)^-(.- 



(59) 



where < > is the mean square local amplitude of 

nearest neighbor spin fluctuations and n is the Bosc func- 
tion. For hole-doped cuprates the g-plateaus constitute a 
problem. For electron-doped cuprates, the plateau width 
is quite small, and the susceptibility is large only over an 
area so the sum rule Eq. 59 is never saturated, and 
Xq and ^ both diverge exponentially with decreasing T. 
For hole doping the plateau width is large, Fig. 4, and 
the curvature on the plateau Aq'^ is relatively small, so as 
T decreases intensity grows all across the plateau. This 
tends to saturate the sum rule, leading to a greatly weak- 
ened divergence of the correlation length, 



a ~ 



ATr\So\ 



UXOQQc 



(61) 



(60) 



(where the latter form follows from Eq. 59, and XOQ = 
XoiQ, 0)) so a ^ at the QCP. 

This result has a number of consequences: (1) neutron 
diffraction^*^'*^ measures the plateau width, Fig. 4 and not 
the correlation length. (2) NMR*^ measures the correla- 
tion length, and in YBCO finds a weak T ^ divergence 
of ^, as predicted. Thus the present results resolve a long- 
standing^^'^^ controversy about the correlation length in 
hole-doped cuprates. (3) In the cuprates, N(;el order ap- 
pears at T > only if the correlation length exceeds 
100a, where a is the lattice constant*^ (the connection 
between Tyv and ^ is discussed in the next section). This 
explains the broad range of hole dopings where there is 
only T=0 AFM order: Fig. 4 shows that the measured 
C 100a only at a: = 0.02 in LSCO (see also Fig. 30, 
below). (4) Moreover, the slope of the r~^/^-term in 
^ decreases rapidly with doping, signalling a QCP just 
above optimal doping. Hence, the triad of features of the 
AFM QCP are also present in the hole-doped cuprates, 
with the broad susceptibility plateaus responsible for the 
striking differences from electron doping. 



C. Incommensurate Magnetism and Competing 
Phases 

The above analysis strongly suggests that at high en- 
ergy scales the physics of the cuprates is dominated by 
magnetic ordering. This includes the large pseudogap 
regime and the attendant QCPs. None of this analysis 
precludes interesting new physics on lower energy scales, 
including of course superconductivity near the QCPs. 
Another possibility is the admixture of a second phase 
generating an enhanced gap - a popular choice being the 
flux phase*^. 

The physics associated with nanoscale phase separa- 
tion, or 'stripe' physics, seems to also fall in this category. 
Incommensurate magnetic modulations are seen in sev- 
eral cuprates - particularly the LSCO family - and while 
the SCR model does flnd an incommensurate suscepti- 
bility particularly for hole doping (Fig. 9b) it probably 
cannot reproduce the observed doping dependence of the 
incommensuration. Indeed, it has been noted^ that the 
incommensurability generally signals an instability to- 
ward phase separation. Experiments suggest that phase 
separation and/or stripe physics is present in the hole 
doped cuprates^^'''^ down to arbitrarily small dopings^^. 
However, the temperature at which stripes are stabilized 
seems too low*^'*'' for them to be directly responsible for 
the pseudogap phenomena. 



19 



A detailed discussion of this issue here is clearly out of 
the question, but the following suggests a possible expla- 
nation. Doped carriers in an AFM are strongly dressed 
by their environment, forming magnetic polarons^" in a 
pure Hubbard model, but in a more general situation 
being sensitive to nearby competing phases^*. Thus, it 
is suggested that the physics of competing phases enters 
the problem at the level of the properties of polarons, and 
different degrees of phase separation and/or stripe forma- 
tion in different cuprates have to do with the tendency 
of polarons to cluster. That is, the stripe physics should 
enter the problem on a lower energy/temperature scale 
than the fundamental pseudogap phenomena discussed in 
the present paper. The g-plateau in hole-doped cuprates 
greatly enhances the sensitivity to stripe physics, since 
the system is close to instability over a wide range of 
incommensurate modulations. 



VII. MAGNETIC PROPERTIES 
A. Electron Doping 

While the present model was developed on the basis 
of ARPES data, the collapse of the Mott gap should be 
clearly reflected in other properties as well, in particular 
in the magnetic response. 

Indeed, Mang, et al.^^ have recently measured the or- 
dered moment M in reduced NCCO samples, and find 
good agreement with the present model^ (see Fig. 30b 
below). The correlation length has not yet been mea- 
sured in reduced (superconducting) NCCO, but there are 
data for the as-grown material, which is insulating®^ 
Fig. 29. The reasons for the striking differences be- 
tween the two types of sample are not fully understood, 
but there seems to be some interstitial oxygen which lo- 
calizes a fraction of the doped electrons, so one must 
dope the as-grown samples more to produce a given re- 
duction of the magnetic properties (e.g., to get a cer- 
tain value of T/v, the doping of the as-grown sample Xg 
must be about 0.02-0.03 larger than for the reduced sam- 
ple Xr, inset in Fig. 29). The data for the undoped 
sample were used to estimate Ps{x = 0,T = 0), and 
thereby u^^ = O.SSAeV. Comparing this to the cr-model 
calculations^^^'^", ps — JS^, gives J — 113meV, in good 
agreement with other estimates. 

However, a fit to Eq. 20 could only be made by reduc- 
ing the (T-dependent) prefactor by a factor of 16. A 
similar problem was encountered in the cr-model calcula- 
tions: one-loop renormalization'^" found ~ SiS 
here (below Eq. 20) , while a two- loop calculation'^^ found 
a T- independent • Introducing a Castro Neto- Hone- like 
interpolation formula^". 
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yields the solid-line fit in Fig. 29, with no adjustment of 
the prefactor. Moreover, the curves for the doped sam- 
ples apply the same correction factors. The agreement in 
T-dependence is quite good; while the theoretical Xr is 
smaller than the experimental Xg, the ratio is consistent 
with both those derived from Tjv and from the magneti- 
zation M, inset in Fig. 29). This strongly suggests that as 
far as magnetic properties are concerned as-grown NCCO 
behaves like reduced NCCO, with a few percent of the 
electrons localized (however, as-grown NCCO never be- 
comes superconducting) . 




FIG. 29. Temperature dependence of correlation length ^ 
in NCCO for X = (solid line), -0.04 (long-dashed line), -0.085 
(dotted line), and -0.10 (short-dashed line). Data are from 
Ref. [ 89]: a; = (open diamonds) and -0.15 (open squares); 
from Ref. [83]: x = (solid diamonds), -0.10 (solid up tri- 
angles), -0.14 (solid down triangles), and -0.18 (solid circles). 
Fits are to Eq. 62, with parameters appropriate to x = (solid 
line), -0.04 (long-dashed line), -0.085 (dotted line), and -0.10 
(short-dashed line). Temperatures are measured in units of 
J — 125meV. Inset: Plot of as-grown nominal doping Xg vs 
reduced nominal doping Xr for fixed values of Tjv (circles) , M 
(squares), and ^ (diamonds). Solid line is Xg = 1.2xr -1-0.012. 

In discussing the as-grown NCCO samples, mention 
should be made of the 'anomalous pseudogap'^^ found 
in an as-grown sample near a:=-0.15 - which should cor- 
respond to a; ~ —0.12 in the reduced samples. From 
Fig. 11c, p, falls off in the range 200 - lOOOK as x 
varies from -0.15 to -0.10, signalling the opening of the 
Mott (pseudo)gap. In the as-grown sample, a pseudo- 
gap was found to open below 240K, centered at 300meV. 
From Figs. 14, 15, the gap near (tt, 0) would be in this 
range. Additional infrared and Raman phonons were ob- 
served, beyond those allowed by tetragonal symmetry. 
This could be associated with the orthorhombic symme- 
try of the magnetic Brillouin zone. Clearly more work 
needs to be done, but if this is the correct interpretation, 
the present model predicts what the doping dependence 
should be, and that similar features should be seen in the 
reduced samples as well. 
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B. Hole Doping 

The results on NCCO should be contrasted to those 
for LSCO^", Fig. 30a, where a saturation of the effec- 
tive ^ is observed in all doped samples. For undoped 
La2Cu04, the data (open circles) largely overlap those of 
Nd2Cu04 (open and filled diamonds), but a small change 
of slope may be present in the best fits. For lightly-doped 
LSCO^o, the data can be fit to Eq. 60 down to ~150K, 
Fig. 30a, below which saturates or decreases. In prin- 
ciple, it should be possible to calculate this saturation of 
^ directly from Eq. 20. As noted in Section III.C, the 
value of A tends to be overestimated when the suscepti- 
bility peak is incommensurate. Thus, the dotted line in 
Fig. 30a is the calculated value of ^, using Eq. 62 with 
parameters appropriate to a:; = 0.10 hole doping, except 
that A/a^ — 0.24, only 1/3 the value estimated from 
Fig. 9a. 

From the a coefficient of Eq. 60 it should be possible 
to extract the magnetization, Eq. 61. However, as ex- 
plained in Section VI, neutron scattering data tend to 
measure the susceptibility plateau width qc , strongly un- 
derestimating ^. Thi s is illust rated in Fig. 30b, where 
magnetization M = V < > derived from ^ via Eq. 61 
is compared to M in NCCO estimated from the ARPES 
data^ (squares) and from magnetization (upright^^ and 
inverted^^ triangles). The ^-derived data include the 
NCCO neutron data of Fig. 30a (triangles) and NMR 
data from YBCO^^ (circles) expected^° to give a bet- 
ter estimate of For both sets of data, the parame- 
ter A/xoqUcl^ was taken as a constant, 2.8. Except for 
the lowest doping, the neutron data lead to an under- 
estimate for M, confirming that the measured ^ is too 
small. In contrast, the NMR data are consistent with the 
electron-doped results, and strongly suggest the presence 
of a QCP just above optimal doping. 
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FIG. 30. (a) Temperature dependence of correlation length 
^ in LSCO for a: = (open circles), 0.02 (inverted trian- 
gles), 0.03 (squares), and 0.04 (triangles) [from Ref. [ 50]], 
compared with Nd2Cu04 (open and filled diamonds, as in 
Fig. 29). Thick solid curve = fit for undoped material from 
Fig. 29). Thin solid lines = fits to Eq. 60. Dotted line = 
calculated value for x — 0.10, as described in the text, (b) 
Comparison of magnetization extracted from £,{T), Eq. 61 
(triangles for LSCO^°, circles for YBCO*^'**"), with that for 
NCCO, taken from ARPES fit, Ref. [ 9] (squares), and from 
magnetization (scaled to M = 0.4 at x=0; triangles: Ref. [ 
92], inverted triangles: Ref. [ 83]. All lines are simply drawn 
to connect the data points, except for that part of the dotted 
line connected with the ARPES data (squares) extrapolated 
beyond x — 0.15. This represents a mean-field calculation, as- 
suming that U does not change with doping over this range, 
and using the band parameters of Ref. [ 9] . 

More recent experiments on very lightly doped LSCO^^ 
have found that the magnetization at these dopings 
is actually incommensurate - consistent with diagonal 
stripes. This points out an interesting parallel with 
the present model: early experimental samples displayed 
flat, diamond-shaped susceptibility plateaus near (7r,7r). 
As sample quality improves, incommensurate structure 
seems to become more prominent: see, e.g.. Fig. 1 of 
Ref. 94. Related behavior arises in the model: The sus- 
ceptibility for hole-doped cuprates displays a flat-topped 
plateau at high temperatures, S > 0. As the tempera- 
ture is lowered, 6—^0, incommensurate structures de- 
velop from fine structure on top of the plateau. Fig. 3 
of Ref. 80, gradually dominating the spectrum. However, 
in the calculations this incommensurability is sensitive to 
sample 'quality': it only shows up when 6 is very close 
to zero. Hence, in real samples, the appearence of such 
structure should be very sensitive to disorder or sample 
inhomogeneity. Finally, it should be noted that inter- 
pretation of the incommensurability in terms of stripes 
remains controversial in cuprates other than the LSCO 
family. Reznik, et al.^^ report an approximately uniform 
ring of incommensurability in optimally doped YBCO, 
which is dispersive and pushed up to finite frequencies 
by the spin gap in the superconducting state. A ring or 
diamond of incommensurability is actually quite close to 
what is found here, and the extension of the present cal- 
culations to the superconducting phase should be quite 
similar to the results of Eschrig and Norman^^. 

VIII. THREE DIMENSIONAL NEEL ORDER 

The (inverse) Stoner factor 5q, Eq. 4, can be general- 
ized to include interlayer coupling: 

Sq{Lu) = S + Aq^ + A,ql - Buj^ -iC'LU, (63) 

leading to a susceptibility 



21 



XQ 



J. 



N) 



f = VT, 



N 



(67) 



(64) 



with ttz = Az/ A. 

In the physical cuprates, the interlayer hopping has an 
anomalous dispersion, generally written as = tzo {cx — 
CyY . This formula holds for bilayer splitting, and in gen- 
eral when the CUO2 planes are stacked uniformly. How- 
ever, as explained in Appendix E, many of the cuprates, 
including NCCO, have a staggered layering, with the Cu 
in one Cu02 plane laying above a vacancy in the neigh- 
boring Cu02 sheet. This leads to a magnetic frustration: 
the Cu in one sheet has four nearest neighbors in the ad- 
jacent sheet, two with spin up, two with spin down. This 
frustration is reflected in a more complicated dispersion 
of tz-. 



tz = tzQiCx -Cyf COS '-^ COS ^ 



a kyG 
cos 
2 2 



(65) 



which vanishes at (tt, 0) and (0, tt), and leads to a greatly 
reduced interlayer coupling. (Effects of AFM frustration 
associated with layering have been discussed in Ref. 97.) 

The consequences of both uniform and staggered stack- 
ing are explored in Appendix E. If the c-axis resistivity 
is coherent, it can be used to estimate the interlayer hop- 
ping tzo- It is found that the value of tzo needed to 
produce a given resistivity anisotropy is approximately 
5 times smaller for uniform stacking, to account for the 
frustration in the staggered stacking. With the corre- 
sponding tzo^s determined from resistivity, both forms 
of interlayer coupling give rise to comparable interlayer 
coupling, and hence a finite Neel temperature. While the 
optimal Q-vcctor depends on doping, at half filling both 
forms predict Q ~ (7r,7r,0), consistent with experiment 
in La2Cu04. Even for quite strong anisotropy, this mech- 
anism can account for the observed T^s (in fact, tends to 
overestimate T/v), without the necessity of invoking ad- 
ditional mechanisms, such as a Kosterlitz-Thouless tran- 
sition, with the reduced spin dimensionality caused by 
spin-orbit coupling effects^^~^°^. 

Within mode coupling theory^"^ (Appendix E), the 
Neel temperature is found from the gap equation 
(Eqs. El, C7) 



Xo{T)U = r,+ 



ttA 



(66) 



where T^d ~ t^ is defined below Eq. E8. It is found that 
T3D is approximately constant, independent of doping in 
the electron-doped regime. Apart from a small numerical 
factor, Eq. 66 differs from the isotropic three-dimensional 
result by the logarithmic factor, which diverges (T/v 0) 
as tz — > 0. 

Equation 66 can be rewritten in a suggestive form. Ap- 
proximating by = A{xoU—ri)/12ua'^ (Eq. 22), then, 
using Eq. 20, the Neel transition occurs when 



where Jz = J(^^o/^)^ J = ^t^/U, and T = AtlQ/UT3D- 
A very similar form was proposed earlier^*^^, and 
experimentally^^ Neel order seems to appear when ^ ~ 
100a. 
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FIG. 31. Comparison of experimental Neel temperatures 
for NCCO and LSCO, (solid line), and for the stripe (mag- 
netic) ordering transitions observed in Nd-substituted LSCO 
[ 46] (solid line with squares) with the model of interlayer 
coupling with staggered stacking and t^o = t/10 ~ 30meV , 
plotted as Tjv/lO (dot-dot-dash line). Also included is the ap- 
proximate expression, Eq. 68 (dotted line with circles). (Note 
that there is a range of hole doping for which A is found to be 
negative; in this range Tjv was arbitrarily assumed to vanish 
in the staggered model, Tjv = 0.) 

Figure 31 compares the calculated value of T/v with 
the experimental values. While the overall doping depen- 
dence is comparable, the calculated Tjv is about an order 
of magnitude higher. The calculation is for staggered 
stacking, with tz adjusted to reproduce the observed re- 
sistivity anisotropy, but Appendix E shows that the over- 
estimate is generic: the coefficient of the logarithm needs 
to be larger to reduce T^. Also shown in Fig. 31 (dotted 
line) is a simplified model, which assumes that 



irA 



3ua^ ln( j|^) 



(68) 



is doping independent, Tq* = 1200if . This model repro- 
duces qualitatively the shape of the numerical calcula- 
tion, but with a magnitude comparable to experiment. 
The magnitude of T/v could be matched almost quanti- 
tatively if Uef / also has a significant temperature depen- 
dence, as discussed in Appendix E . The overall doping 
dependence is also comparable to experiment. The agree- 
ment could be further improved by using a smaller value 
of t' , which would shrink the doping range over which 
Neel order occurs. 

Finally, it should be noted that a finite T/v can change 
a continuous QCP into a, first order. This follows beacuse 
the plateau width increases with increasing temperature. 
Hence, near the plateau edge, the system can satisfy the 
Stoner criterion at some finite temperature, but fail to 
satisfy it at a lower temperature, having fallen off of the 
plateau edge. Such a first order termination of the AFM 
state seems to be found in the electron doped cuprates. 
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most notably in Pr2_a;Cea;Cu04 (PCCO)"'^'''*, and in a 
related organic material^°^. 

IX. DISCUSSION 
A. Slater vs Mott Physics 

Theories of magnetism fall into two diametrically op- 
posed classes^°^: band vs atomic models, or Slater vs 
Mott physics. In both approaches, the copper band is 
split by a Mott gap into two parts. In Slater theory, long 
range magnetic order leads to a unit cell doubling, so 
each subband remains conventional, with two electrons 
per unit cell. In Mott theory, the bands are highly un- 
conventional: the gap opening is purely a local effect - 
there is an energy penalty of U for two electrons to sit 
on the same copper site. Since there is no change in lat- 
tice symmetry, the unit cell remains the same, and the 
bands hold only half as many electrons as conventional 
bands. Distinguishing between the two models is com- 
plicated, since in the Hubbard model, residual hopping 
proportional to At'^/U ^ J leads to AFM coupling of the 
electrons, and can lead to parasitic Neel order, at a tem- 
perature T/v much lower than that at which the Mott 
gap opens. On the other hand, strong fluctuations in the 
Slater model can greatly reduce Tn, leaving a pseudogap 
near the mean-field instability temperature. 

While the mean field results might give a good quali- 
tative picture of the Slater regime, they are unlikely to 
be able to describe Mott physics. However, it would still 
be hoped that the mean field results can give an indica- 
tion of when the crossover is likely to occur. Here, two 
separate indicators are presented. 

A first indication comes from looking at competing or- 
ders. A Stoner criterion Uxq — ^ gives the onset tem- 
perature for magnetic order at q, ranging from AFM, 
q = Q to ferromagnetic, q = 0, Fig. 32. While at half 
filling for any value of U, AFM order dominates, the 
splitting decreases with increasing U. The local, or Mott 
physics should arise when fluctuations to all magnetic 
orders are comparably likely, or the spread in transition 
temperatures, ATc, is << T. Since the probability of a 
fluctuation of N particles into a phase with excess free 
energy A/ is ^ ^-NAf/ksT^ ^^^^ ^^^^ crudely state that a 

phase will be significantly excited if T^v — Tc{Q) < aoT/v, 
where ag is a small numerical constant. The width of the 
ATc curve in Fig. 32a shows the fraction of the Brillouin 
zone that is significantly excited for ao = 0.01 (e.g., for 
U/t > 32.5, all modes are excited). This suggests that 
for U > 15i, these fiuctuations spread over a significant 
fraction of the Brillouin zone, while for U > 30t virtu- 
ally all magnetic states are equally excited and the Slater 
picture is badly broken down. However, the cuprates are 
generally found to be in the regime U < 12t, where a 
Slater picture should be reasonably accurate even close 
to the Ttr crossover. 
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FIG. 32. (a) Mean field transition temperatures for Neel 
(Tjv) and ferromagnetic (Tc) orders and their difference 
ATc = Tn — Tc (upper dashed line). At any U/t the ra- 
tio of the shaded area to the total area below this line gives 
the fraction of the Brillouin zone which is significantly excited 
{Tc{Q) > (l-ao)Tiv, for ao = 0.01). (b) Replot of transition 
temperatures scaled to U. (c) Plot of Tc vs q, for U = 6t. The 
curve bears an uncanny resemblance to the (scaled) electronic 
dispersion of the LHB, long-dashed line. 
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FIG. 33. (a) Bare susceptibility xo(Q,0) at T = IK, for 
several values of r = -0.0552, -0.110, -0.276, -0.552, -0.828, 
-.9, -.99 (solid lines), and (dashed line), (b) Crossover cou- 
plings as a function of r: Uv (triangles), U2 (squares), and 
Ui (circles). (Dashed line — Ui for electron doping.) 

Alternatively, when the mean-field solution becomes 
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insensitive to the band structure, it is likely that a local 
picture is becoming dominant. In the present instance, 
the band structure is determined by the ratio t' /t. For 
any non-zero value of t' , the susceptibility has a generic 
doping dependence^°'''^°, Fig. 33 - changing the sign of 
t' merely interchanges electron and hole doping. The 
role of the susceptibility plateaus can be quantified, by 
defining ranges of U where the nature of the transition 
changes, Fig. 33b. Thus, for U < Uy, there is no Mott 
transition at x = 0, and the physics is dominated by an 
AFM transition at the VHS; for Uy < U < U2, there is a 
Mott transition at half filling, which terminates (on the 
electron-doped side) before the plateau ends, and hence is 
controlled by dynamic critical exponent z = 2; for U2 < 
U < Ui, the Mott gap collapses in the enhanced regime 
near the edge of the plateau; and for U > Ui the Mott 
gap terminates well off of the plateau, in a region of z ^ 1 
physics. For the present r = —0.552, the approximate 
values are Uy/t = 2.4, U2/t = 2.6, and Ui/t = 3.6 {x < 0) 
or 4.1 {x > 0). Note that the z — 2 regime is quite 
narrow, and can probably be subsumed into the VHS 
regime. These values depend on t' , and the VHS moves to 
half filling as t' — > 0, Fig. 33b. Even when t' = 0, Sen and 
Singh^°* find a crossover from SDW-like to Heisenberg- 
like behavior as correlations increase beyond Uov = 3.26t 
(diamond in Fig. 33b). It must be kept in mind that U 
depends on doping, and the above estimates refer to U 
near the plateau edge. The bare Uq = U[x — 0) can 
be estimated by assuming the doping is high enough to 
reach the Kanamori limiti'\ U = C/o/(l + Uo/8t). This 
results in Uov/t = 3.4, C/02A = 3.8, Uoi/t = 6.5 or 8.4 
for T = —0.552. These last values are comparable to 
but somewhat smaller that those estimated by the first 
criterion. 

Note that the cuprates are in the range U2 < U < Ui, 
where the plateau edges form natural phase boundaries 
for the Stoner criterion, thereby providing a natural ex- 
planation for the approximate electron-hole symmetry of 
the QCPs. 

From the above discussion, it might be concluded that 
the cuprates are in the Slater regime, where mean field 
results are qualitatively accurate. While this is probably 
true for the electron-doped cuprates, it must be kept in 
mind that the criterion for the breakdown of the Slater 
regime is the flatness of the susceptibility in q - the inabil- 
ity of Slater theory to deal with many competing phases. 
In this light it is interesting to speculate whether the flat- 
ness of the g-plateaus might be a signal of enhanced Mott 
physics in the hole-doped cuprates. 

B. Magnon Bose Condensation and Non- Fermi 
Liquid Physics 

Figure 34 shows the sharp peak which arises in /mE at 
low T. The growth is exponential, approximately match- 
ing that of the coherence length, Eq. 20. (Note that it 



requires a fine mesh in the integral of Eq. 28 to capture 
this growth.) This peak arises exactly at the incipient 
magnetic zone boundary, and turns into true Bragg scat- 
tering at the transition to long range order: the increase 
in peak height is almost exactly compensated by a de- 
crease in the width of the peak. A simple physical expla- 
nation is that the SDW transition can be interpreted as a 
Bose condensation of the zone boundary magnons. Then 
the Mermin- Wagner theorem reduces to the fact that in 
a two-dimensional system, Bose particles can only con- 
dense at r = 0. A similar explanation for the transition 
has been presented earlier'^. 
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FIG. 34. (a,b): Blowups of /m(E) for a; = at (7r,0) (a) 
and (7r/2,7r/2) (b) at T = lOOK (solid lines), 500K (long 
dashed lines), and lOOOK (a) or 750K (b) (short dashed 
lines), (c) Maximum of /m(S) vs T for (vr, 0) (squares) 
and (7r/2,7r/2) (circles); 0.1/(full width at half maximum) 
for (7r,0) (triangles) and (n/2,n/2) (diamonds); solid line — 
corresponding ^(T), Eq. 20. 

In turn, the soft zone-boundary phonons explain one 
origin of non-Fermi liquid physics in the model: Bragg 
scattering from a fluctuating diffraction grating. How 
does one deflne Luttinger's theorem when the unit cell is 
strongly fluctuating? 

C. Comparison with Other Calculations 

As noted above, the present calculations predict that 
at the magnetic QCP the (possibly T = 0) Neel phase 
will terminate, the Mott gap will collapse, and the Fermi 
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pockets will merge into a large Fermi surface. This result 
distinguishes the present calculations from many others 
in the literature. Here a number of Slater-like theories 
are discussed. 

The present results are generally consistent with the 
t — J model in the low doping regime. However, since 
the t — J model cannot readily deal with both Hubbard 
bands simultaneously, it is not appropriate in the present 
analysis of electron doped cuprates. where (a) ARPES 
can detect both bands (at least up to the Fermi level), 
and (b) the Mott gap is found to collapse with doping, 
leading to an overlapping of both bands at the Fermi 
level. In the t — J model double occupancy is forbidden 
in the LHB, while in the UHB empty sites are forbidden. 
Moreover, the Hubbard model only allows J values for 
J < 1 (Figs. 25, 27), and near this upper Hmit signifi- 
cant modifications are needed. In the SCBA approach to 
the t — J model, the parameter A\ has a broad peak^^ 
when J ~ 0.8, not found in the Hubbard model SCBA 
calculations, suggesting that t — J and Hubbard can be 
equivalent only for J << 0.8 or for {/ >> 5 - that is, 
near half filling only. It is interesting to note that a re- 
cent t — t' — t" — J model calculation seems consistent 
with the first doped carriers forming weakly interacting 
quasiparticles in pockets of the respective upper or lower 
Hubbard bands, for either electron or hole doping^°^. 

The NAFL and spin fermion models are also based on 
Slater- type physics, and should in principle make similar 
predictions to the present SCR model. However, they 
tend to take their parameters from experiment, which 
can lead to complications in the presence of stripe phases. 
For example, for hole doping, long range Neel order and 
diverging susceptibilities terminate at a very low doping, 
X ~ 0.02. While the SCR model predicts xqcp ~ 0.25, 
some empirical models take xqcp ~ 0.02. In this case, 
the QCP is divorced from Mott gap, since the Mott gap 
will clearly persist above x = 0.02. Even worse. Matsuda, 
et al.*® have shown that for doping between x — 0.02 
and half filling the system is phase separated, so uniform 
AFM order exists only at a; < 0. 

Three examples of spin fermion calculations will be 
given, to highlight the differences and similarities. (1) 
Abanov et al.^^° postulate a small-x magnetic QCP. They 
find that the magnetic resonance mode frequency goes to 
zero at this QCP, but also the superconducting gap van- 
ishes at the same doping, which would have important 
consequences for the mechanism of superconductivity. 
This is in sharp contrast to the present model, where the 
magnetic QCP is at much higher doping. In this model, 
the I0W-.T QCP is superconducting, presumably associ- 
ated with stripe effects. (2) Chubukov and Morr^^'^^^ 
studied the crossover from small to large Fermi surfaces - 
here driven at fixed doping by reducing U. They state^^^ 
that, "as the system moves away from half filling, the 
spectral weight transfers from the upper band into the 
lower band and, near optimal doping, there exists just 
one coherent band of quasiparticles." This suggests that 
the crossover is due to a spectral weight shift, and not to 



the gap closing - in contrast to the present results and 
to experiment on NCCO. However, it should be noted 
that (a) their paper actually concentrates on changes at 
the Fermi level, and did not explore how the UHB might 
have shifted with U; and (b) it is possible that the co- 
herent part of the UHB collapses, while some weight re- 
mains in the incoherent part. (3) On the other hand, 
Schmalian, et al.''^' go beyond the SCBA, summing both 
non-crossing and crossing diagrams via a generalization 
of a technique of Sadovskii^^^; their results for hole dop- 
ing are quite similar to the present results, with a mag- 
netic QCP above optimal doping but with ^ adjusted 
at each doping to fit the experiment. 

A number of groups have studied the Hubbard model 
using FLEX calculations, and have had considerable 
success in describing anomalous transport properties^^^. 
Here a pseudogap is found even though the FLEX model 
cannot describe the splitting into UHB and LHB, and 
the pseudogap is derived from superconducting fluctu- 
ations. However, these models are consistent with the 
present results, in that (1) the pseudogap they describe 
is clearly the lower, leading edge pseudogap which is not 
described by the present model, and (2) their calcula- 
tion of the normal state properties require a value of 
U/t ^ 1.5 — 2.5 much smaller than the values found at 
half filling, and comparable to (or even smaller than) the 
doped values found here. [Spin fermion calculations also 
extract a small value of U there called g - from exper- 
iments in near-optimally hole doped cuprates^^*.] 

The present calculations are in general consistent with 
the results of Rcf 3. These authors employ a (Two- 
Particle Self Consistent) conserving approximation, and 
attempt to calculate U{x) directly^^^. However, they in- 
corporate the strong thermal (Mermin- Wagner) fluctua- 
tions directly into their deflnition of [/, so the resulting 
doping dependence should not be compared to the form 
assumed here. 

A leading edge pseudogap can also arise in the Hub- 
bard model in the absence of superconductivity^^^'^'''^^, 
but only for large U > 8t^'^. 

Some recent calculations have conflrmed that U must 
decrease with electron doping to reproduce the ARPES 
data: in Kusonose and Rice^'^ the demonstration is indi- 
rect - the gap collapse does not occur in a SCBA calcu- 
lation if U is kept large. Scncchal and Tremblay^** give a 
more direct demonstration; their model can also explain 
the hole-doped pseudogap near (7r,0) in the absence of 
stripe physics if U does not decrease with hole doping. 

Finally, a proper study of the model incorporating 
QCP fluctuations is a strong desideratum, but the prob- 
lem of combining QCP and Mermin- Wagner fluctuations 
has rarely^^'' been tackled in the literature. 
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D. VHS 

Whether or not the VHS is responsible for the observed 

electron-hole asymmetry, the present calculations reveal 
some novel features of Van Hove physics. 

1. Temperature Dependent VHS 

As noted by Onufrieva and Pfcuty^^ , the VHSs associ- 
ated with the susceptibilities (and hence with charge or 
spin nesting) are different from those associated with the 
density of states (and superconductivity). Thus, whereas 
superconductivity will occur at the same optimal doping 
for all temperatures, the doping of maximal nesting in- 
stability is a strong function of temperature. 

This contrasting behavior of nesting vs pairing suscep- 
tibilities is related to a characteristic difference in the 
nature of the two instabilities. A superconducting in- 
stability has an intrinsic electron-hole symmetry, which 
means that the gap is tied to the Fermi level, and a full 
(s- or d- wave) gap can be opened at any doping level. 
On the other hand, a nesting gap is dispersive, and only 
part of it lies at the Fermi level (except in special cases). 
Furthermore, a (superlattice) Luttinger's theorem must 
be obeyed, requiring the presence of residual Fermi sur- 
face pockets. Stated differently, a full nesting gap can 
only open at integer filling, so as the interaction strength 
increases, any nesting instability must migrate to inte- 
gral doping (e.g., half filling in the original band struc- 
ture). This same VHS migration is mirrored in the T- 
dependence of the magnetic (or charge) susceptibility. 

2. VHS Transitions 

We have seen that the doping-dependent Ueff gives 
rise to a Mott gap collapse near the edges of the suscep- 
tibility plateau in Fig. 1. If Ueff is smaller (dot-dashed 
line: Ueff reduced by 2/3), more complicated behavior 
should arise. Due to the peak in x near the H-point, 
there could be a reentrant transition, with one magnetic 
order near half filling, and a second near the VHS. For an 
even smaller Ueff (or replacing Ueff J)'^^, the transi- 
tion near x = can be eliminated, leaving a spin density 
wave transition near the VHS. In principle there could 
even be a phase separation between two AFM phases: an 
insulating phase near half filling and a metallic phase 
near the VHS. 



X. CONCLUSIONS 

The key conclusion to this work can be stated as fol- 
lows: In doped cuprates there is a magnetic QCP where 
three factors coincide: the crossover from small to large 
Fermi surface, Mott gap collapse, and Neel transition 



termination. In the SCR calculation there is no finite 
temperature Neel transition, at least in the isotropic 2d 
limit, but the zero-temperature Neel transition persists 
with doping up to a QCP controlled by a modified Stoner 
criterion. While the Mott gap opening is more of a 
crossover than a sharp transition, nevertheless, the up- 
per and lower Hubbard bands merge at nearly the same 
point, and the Fermi surface pockets recombine to a sin- 
gle large Fermi surface, consistent with band structure 
calculations. Comparison with experiment suggests that 
this correctly describes the situation in electron doped 
NCCO, both in ARPES (Section IV) and in magnetiza- 
tion studies, Ref. 83 and Section VII. 

The hole doped case also appears to fit this model, but 
with complications associated with the g-plateau. Thus, 
(a) the pseudogap collapses in a QCP, as expected; (b) 
evidence for the Fermi surface crossovcir has recently been 
reported^^'^*; (c) the correlation length appears to di- 
verge as T ^ 0, but much more weakly than for electron 
doping, due to a sum rule saturation***^. 

In more detail, the main results of this paper can be 
summarized: 

• Fluctuation effects were added to the mean field Hub- 
bard model via a mode coupling calculation, which al- 
lowed satisfying of the Mermin- Wagner theorem {Tn — 
0). It was found that the mean- field gap A^f and Neel 
temperature evolved into a pseudogap Ap^ ~ ^m/ 
and an onset temperature T* ~ T™-'^ (as is familiar from 
the related CDW resuhs). 

• The resulting dispersions and Fermi surfaces are in 
excellent agreement with photoemission experiments on 

electron-doped cuprates^, while the pseudogap seems 
consistent with ARPES and tunneling results in hole 
doped cuprates. 

• Magnetic properties - saturation magnetization and 
coherence length - are also well fit by the same model. 
The good agreement between ARPES and direct mag- 
netic measurements leaves little doubt that the (large) 
pseudogap is predominantly magnetic in origin. 

• The zero-temperatme Neel transition is controlled by 
a Stoner-like criterion, hence is sensitive to the bare sus- 
ceptibility and in turn to the Fermi surface geometry (hot 
spots). This lead to an approximately electron- hole sym- 
metric QCP near optimal doping (termination of hot spot 
regime), at which both zero temperature Neel transition 
and pseudogap transition simultaneously terminate. 

• The model leads to a NAFL-type susceptibility, and 
the calculation of the NAFL parameters has been re- 
duced to a calculation of the coupling parameters U and 
u, the former having a significant doping (and possibly 
temperature) dependence. At present, U{x) is estimated 
from experiment, and the mode coupling u via consis- 
tency with the t — J model. (A small portion of the 
renormalization of U arises from quantum corrections to 
the Stoner criterion.) 

• The present theory differs from conventional NAFL 
theory by the inclusion of two cutofi^ parameters, Qc and 
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uj~ , which shrink to zero at either the H- or C-points. 
For example, qc is large near the H-point, but shrinks 
to zero at the C-point, causing the A parameter to have 
a strong temperature dependence in the electron-doping 
regime. 

• Finally, a striking temperature/frequency dependence of 
the VHS susceptibility peak'*^, causing it to shift to half 
filling at high T, is interpreted in terms of Luttinger's 
theorem: if the coupling is strong enough to open a full 
gap, the gap must fall at half filling. 

Note: After the present work was completed, I re- 
ceived a preprint from A. -M.S. Tremblay reporting simi- 
lar calculations for electron-doped cuprates^^^. 
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APPENDIX A: THREE BAND MODEL 

A major simplification of the present calculation is to 
treat the cuprates in a one-band model. This is consis- 
tent with the Zhang-Rice picture-'^^^, although the ap- 
proximation is less drastic for electron doping, since the 
upper Hubbard band is already predominantly copper- 
like. Nevertheless, the model also describes the doping 
dependence of the 'lower Hubbard band', which is really 
a charge transfer, predominantly oxygen-like band. Here 
an explanation for why this simplification works is sug- 
gested. 

Even without carrying out self-consistent calculations, 
the nature of the Mott transition can be understood by 
introducing a doping dependent gap. The energy bands 
can be calculated from the hamiltonian matrix 

H = J2 ^d]dj + ^ tcuo[d]p^ + (c.c.)] 

<ij> 

+ too[p]pf + {c.c.)] 

<3J'> 

+Unjinji + UpUi^nn), (Al) 

where A is the difference in on-site energy between cop- 
per and oxygen, tcuO is the copper-oxygen hopping pa- 
rameter, too the oxygen-oxygen hopping parameter and 



U (Up) the Hubbard interaction parameter on Cu (O). 
For good agreement with the doping dependence of the 
one band model, it is necessary to properly incorporate 
the Hartree correction to the self energy, A = Ao -I- S/f , 
Sjj = Un^ (for up spins), and = n/2 — mg, with n 
the average electron energy. The resulting dispersions are 
shown in Fig. 16 for the antibonding bands, and Fig. 35 
for the full dispersion. In these figures, the following pa- 
rameters are assumed: icuO = 0.8eT^, too = —Q.AeV, 
Ao = 0, J7 = 6ey, and Up = i.lbeV . 

The band dispersion is extremely similar to that 
found in the one band model. Fig. 15, even though the 
lower band crosses over from the Zhang-Rice (hybridized 
copper-oxygen band) at half filling to a more copper like 
lower Hubbard band with increasing electron doping. In 
addition, the effective magnetizations are proportional. 
Fig. 36, although the one-band model overestimates the 
magnetization by 1/3. This can be understood: in the 
three-band model, the shape of the Hubbard bands is 
fixed by the combined effects of the magnetic instability 
and hybridization with the oxygen band. In the one band 
model, only the former effect is present, necessitating a 
larger value of m to produce the same net splitting. 
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FIG. 35. Dispersion of all six bands in three-band model, 
assuming mq = 0.3 (a) and 0.01 (b). 

This remarkable agreement between one and three- 
band models goes well beyond the Zhang-Rice model. 
That model is restricted to the LHB in a small range of 
doping near half filling; the present results compare both 
LHB and UHB over the full range of electron doping. 
The result is nontrivial - in the three band model, the 
bonding and non-bonding bands are also split into up- 
per and lower Hubbard bands. This degree of agreement 
comes about because the parameter A includes a large 
contribution from the magnetic Hartree term. In turn, 
this suggests that in the absence of magnetic effects the 
Cu and O energies are nearly degenerate - as found in 
early LDA band structure calculations (see discussion in 
Ref. 120). 
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FIG. 36. EfTective magnetization me// = mU/Qt for the 
three-band (circles) and one-band (squares) models. The 
one-band result has been multiplied by 3/4 to better agree 
with the three-band results. 



APPENDIX B: CHARGE SUSCEPTIBILITY AND 

Ueff 

The proper choice of vertex corrections is an unre- 
solved issue in the analysis of the Hubbard model. It 
is known to be of critical importance for generating a 
pseudogap^^^. Here, by comparing simple mean- field and 
SCR models to experiment, it is shown that the net ef- 
fect of vertex corrections is to make the coupling U ef- 
fectively doping (and possibly temperature) dependent. 
Kanamori^^ showed that the effective Hubbard U should 
decrease with doping, as an electron can hop around, 
and hence avoid, a second electron. In the limit of a 
nearly empty (or full) band, this should lead to a correc- 
tion of the form Ueff ^ U/{1 + U/W), where W = U 
is the bandwidth. It was found^^^'^'^ that Monte Carlo 
calculations of the susceptibility of a doped Mott insula- 
tor were approximately equal to the RPA susceptibility 
with suitable ?7e//, and Chen, et al.^^^ suggested the ex- 
plicit form Ueff = U/{1+ < P > U), with P given by 
a vertex correction to the susceptibility and < • • • > an 
average over q, at zero frequency. Figure 37b presents 
a calculation for Ueff based on Chen, et al. However, 
whereas Chen, et al. performed the average in the para- 
magnetic phase, using bare Green's functions, here the 
dressed Green's functions appropriate to the Neel phase 
are used, to approximately incorporate the effect of this 
gap. This makes little difference, since P is dominated 
by the intraband terms, and remains finite at half filling. 
Explicitly, 



^-fl-fl' 



E^{k) + Ej{k') - u) - i5 



E±{k) = -{ek + ek+q±En), 



(Bl) 

(B2) 
(B3) 



U^^,{k,k') = -{l + iAk){l+jA,,,) + ijBkBk', (B5) 

with i,j summed over +, — , A the AFM gap, and Ak — 
{ek-ek+Q)/Eok, Bk = A/Eok- In agreement with Chen, 
et al., the calculation finds U to be renormalized by a 
factor of 2 at finite doping, but does not recover a large 
U near half filling, althovigh different results are found 
depending on whether a; = from the start (triangle) or 
whether x — > from the hole or electron doping sides. 




0.1 0.2 0.3 0.4 



(b) 



0.1 0.2 0.3 0.4 



FIG. 37. Calculated Ueff assuming (a) simple screening or 
(b) full vertex correction of Chen, et al. [ 122]. In both cases, 
a bare U — 6.75t was assumed. Solid lines — electron doping; 
long dashed lines = hole doping; triangles (squares) in a = 
paramagnetic screening of f/, at T = IK (2000K); triangle in 
b = undoped; circles = data of Ref. [ 8] . 

For modelling purposes, it is useful to have a Ueff 
which evolves smoothly from a large value at half filling to 
a reduced, Kanemori value at finite doping. A simple toy 
model consists of taking the RPA screening of a charge 
response. There should be a close connection between the 
Kanemori mechanism and screening. Screening involves 
creation of a correlation hole about a given charge, while 
Kanemori's Ueff involves the ability of a second charge to 
move around the first, while avoiding double occupancy. 
Near half filling, the second charge must move in the 
correlation hole. Approximating^ the vertex correction 
by the RPA screening of the charge susceptibility. 



Ueff = 



U 



i+<x>u' 



(B6) 



it is possible to reproduce^ the experimentally observed* 
doping dependence, while matching the calculation of 
Chen, et al. away from half filling. Fig. 37a. 

In this calculation, issues of self-consistency are also 
important. To minimize screening at half filling, it is nec- 
essary to reproduce the gap in the susceptibility. Hence, 
the susceptibility in Eq. B6 is approximated by the charge 
susceptibility in the AFM state, Xo° from Eq. 2.24 of 
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Rcf. 1, evaluated with the bare U = 6.75t. (In princi- 
ple, at finite doping there is a coupling to the longitu- 
dinal magnetic susceptibility^'^, but this is neglected for 
simplicity.) The importance of using the AFM suscep- 
tibility is illustrated in Fig. 37a: the solid and dashed 
lines show Uejf calculated using the charge susceptibil- 
ity in the Neel state, while the corresponding lines with 
triangles use the paramagnetic susceptibility at low T. 
The latter calculation finds a nearly doping independent, 
but small Ueff, the former reproduces a large, weakly 
screened U near half filling. Such a difference is expected 
in terms of screening: when there is no gap at half fill- 
ing, the enhanced susceptibility should be better able to 
screen U, resulting in a smaller Ueff - This suggests that 
Ueff should have an important temperature dependence 
as the gap decreases - which in turn will cause the gap to 
close at a lower temperatiirc. Figure 37a also shows that 
there is a weak temperature dependence of the screening. 
The calculations suggest that the large values of U found 
in the cuprates are characteristic mainly of the half filled 
regime and relatively low temperatures. A similar but 
larger screening effect was recently reported by Esirgen, 
et al.i23. 

This procedure is still not fully self consistent. If there 
is a large difference between the bare U and the screened 
Ueff, the gap in x should depend on the actual Ueff. 
However, since Ueff ^ U at half filling, any simple im- 
provement will not significantly change the overall dop- 
ing dependence. This is the same kind of lack of self- 
consistency found for the SCR approach, and will be here 
neglected. 



APPENDIX C: IMPROVED SOLUTION OF SCR 
EQUATION 

Approximating coth(x) = max[\lx,\)^ and introduc- 
ing the notation Aq^ = Aq^ +5, Uq = Aq'^/a^^, and 
t = 2TC, the solution to Eq. 6 becomes 



with 



= *[^i(^)-^i(|)]' (C3) 



h{x) = Io{tan~^{x)) - tan~'^{x)ln{x), (C4) 



In{tan0)d9 = L{x) + L{- - x) - L(-), 

(C5) 

= — Jq ln{cost)dt is the Lobachevskiy 



and L{x) = 
function^^^. 

For most purposes, it can be assumed that 6 «t « 

Aq'^, fto), in which case Io{tan~^{x)) = 9{ln (6) — 1), with 
9 = min{x, 1/a;}, and then F2, Eq. C3, simplifies. 

F2 = ln{*-)[S + ttan-\*-)] +^-^ 

^^tlnii). (C6) 
Defining Z =1 + {3ua'^ /ir'^AC)ln{a^/t), then 

Z5-So = 7-ln{^), C7 

nA 

which agrees with Eq. 13 when Z 1. 



APPENDIX D: PARAMETER EVALUATIONS 

At T = 0, the imaginary part of the susceptibility 
x((5, w) can be calculated analytically: 

/m(x(<3, w)) = Y,^f{ek) - fi^k+Q))^i^UQ -^k-^) 

k 

F{9i,k) - Fi92,k) 



At 



(Dl) 



3ua^ 



with 



S+Aq^, f'O^ 

dy / dx 
S Jt 

Aq^ 



x^ + 



'■-ln\ 



0? + {t/a^Y 



a^tan ^(a,) — -ln[ 



5, + 



F2=t dy 



a a;2 + y2 



f'+^'i'^dy ,,t, 
= t / — ian"V-) = 
Js y y 



(CI) 



2^..,,^^J-a.tan-(-) (C2) 



where F{6,x) is an elliptic integral, k = s/l — {u/8t)^, 
and sin{6i) = sin{(j)i)/k, with 



COS^ ((/>l) 



if a; < 
uj/2 if a; > I 



2 /- J, \ J if w < Wn 
cos^ (^2) = i 1+ -t Z ' 
[1 if w > Wo 

with /t = /i/2f, Cj = (jj / At, Cj. = a± - 



(D2) 



(D3) 



and 



a± = 1 — (/'},±c2')/r. Similar results for t' = are discussed 
in Rcf. 40. The real part Rex can be found from the 
Kramers- Kronig result. 



1 r°° 

Rex{Q,oj) = - / 



Imx{Q,i^')oj' dw' 



(D4) 
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1. A at the C-point 

To understand the q-plateau, and in particular the C- 
point, where the plateau width shrinks to zero, it is con- 
venient to introduce a simplified model^^^, for which the 
g-dependence of x can be calculated analytically. While 
the dashed lines in Fig. 41 represented an uj shift, they 
can equally well describe the g-shift of the energy denom- 
inator, Eq. 3. The plateau edge corresponds to the point 
where the dashed line intersects the q-shifted FS (hori- 
zontal arrows). (Recall that q = Q + cf .) In the simplified 
model, the energy denominator is linearized, so Ae oc k±, 
independent of fc|| . Chosing q to point along the (tt, tt) 
direction, the FS can be approximated by two circles of 
radius kp, centered at (tt, tt) and (— tt, — tt) (for this choice 
of q the other two circles at (tt, — tt) and (— 7r,7r) can be 
ignored). The Q-shifted FS is then a circle centered at 
F = (0,0). The FS at (tt, tt) and the Q-shifted FS are 
illustrated in Fig. 38c. To keep the picture symmetri- 
cal, both FSs are shifted (in opposite directions) by q'/2 
when q' 0. 
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FIG. 38. (a) Calculated susceptibility xil) for several val- 
ues of overlap 5. (b) Blowup of plateau region, for Xg ~Xq=o- 
(c) Model of Fermi surfaces, defining 5, ki (fcj_) and 

Adding the contributions of the overlap of the q-shifted 
FS with both the FS at (tt, tt) and the one at (— vr, — vr), 

Xq oc Iks+q'/2 + Iks-q'/2, with 



'^^ dk 1 dk 



(D5) 



where the region of integration is over the part of the 
upper FS in Fig. 38c not overlapped by the lower (q- 
shifted) FS, and k± ranges from zero at the apex of the 
wedge to the middle of the upper FS, kc = kp — ks, where 
ks is the overlap parameter defined in Fig. 38c. To lowest 
order, for kg << kp, 



I = 2kp + \/kpks In 



1-/3 , 
1 + 



(D6) 



with f3 = ^/ks/kp. The expression for Iks-q'/2 must 
be modified when q' > 2ks and the two FSs no 



hs-q'/2 



2kp[l - 7tan-i I/7], with 
7 = a/ {q' — 2ks)/2kp. The calculated susceptibilities, 
Fig. 38a, display the flat topped plateaus with weak pos- 
itive curvature {A < 0, Fig. 38b). At the plateau edge 
the susceptibility falls sharply, x~ 1 — 7r7/2~ y/q'. The 
C-point corresponds to fc^ = 0. 



2. C and Plateaus in Frequency 
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FIG. 39. (a) Imx{Q,uj), (b) RexiQ,uj), (c) Imx/u; = C, 
and (d) dRexiQ , ^) / duj , for (s-ic): H = (solid line), -0.05 
(long dashed line), -0.10 (dashed line), -0.15 (dotted line), 
-0.20 (dot-dashed line), -0.25 (dot-dot-dashed line), and 
-0.30eV (short dashed line); (b,d): a; = (solid hne), 0.04 
(long dashed line), 0.10 (dashed line), and 0.15 (dotted line). 

Figure 39 illustrates Imx{Q,i^), Rex{Q,uj), and 
Imx/i^ = C. While plateaus in i?e(xo) have been 
noted above. Fig. Ic, here the main interest lies in C = 
Ulmx/i^- This linear-in-frequency contribution to Imx, 
generated by hot spots, is an important parameter in 
SCR and NAFL theories, and has been well studied. The 
height of the plateau at zero frequency C = UC{uo — 0) 
can be represented as a frequency uji — I/C, which can 
be found explicitly'^^ 



C = 



2ttJsIq{1 + rc^o) 



1 



(D7) 



(with J = U^/U,_ sl, = l- cl,). 

However, it is important to note that C also approxi- 
mates a plateau, particularly near the H-point, Fig. 40, 
with a well-defined cutoff. Moreover, the width of this 
plateau vanishes near both the H- and C-points, con- 
trolled by two characteristic frequencies, w^T, Eq. 24, and 



8t, 



^0 = — [v^l - At - 1], 

T 



(D8) 
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respectively. The origin of these critical frequencies can 
be understood from Fig. 41. The thick (thin) sohd hues 
represent the original (Q-shifted) Fermi surfaces, while 
the dashed lines represent 



for various values of lo. Equation D9 gives the points at 
which the denominator of Xo(Qi Eq. 3, vanishes. Thus 
at T = 0, Im{xo{Qi^)) is proportional to the length of 
the dashed line lying between the original and Q-shifted 
FSs (i.e., where f {ej:) - f {^k+q) = ±1)- Since the two FSs 
meet at an angle, forming a wedge, /m(xo(Q, t<^)) ~ ^■ 
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FIG. 40. C calculated for several values of /i = -0.355 
(diamonds), -0.357 (circles), -0.358 (squares), -0.359eV (trian- 
gles) = -0.3599eV]. Inset: Band dispersion eg (solid line) 
^k+Q (dashed line) , for /i = 0. Arrow = uj~ . 



set. Fig. 40), while near the C-point it is wq- The ver- 
tical arrows in Fig. 41 indicate cj^, where the dashed 
line (Eq. D9) intersects the FS at the zone boundary, 
while the horizontal arrows are**' loq, where the dashed 
line ceases to intersect the Q-shifted FS. There is a 
crossover at /ic — — 0.14ey: for ^ > fic, — lo~ 
while for /i < /ic, lo^) — lo~ . Combining Eqs. 24, D8, 
ujQ = uj- at [ic = [l-z(2- Vi)2]2t/T = -0.1384ey, with 
z — \ — T. For cj > mm{cJo,Wj7}, /m(xo(Q,t^)) ^ cj-^/^, 
so C ~ 1/w^/^ - i.e., the susceptibility is no longer on 
the plateau. 

Defining a width parameter ol^jj = mm{a~,a2,}, with 
Q!° = wo/t^i, then 

_ 27rt(l-r) 1 + rc^On 1 .-p,.„^ 
^i/^c = 77 [ ] = — • (DIO) 

U -r 

This latter is in good agreement with the numerical re- 
sults (arrows in Fig. 41) and is similar to the result found 
by Onufrieva and Pfeuty'*^, using a hyperbolic band ap- 
proximation valid near a VHS, wi/w^T — 27ri(l — t)/U. 

Because of the dynamic scaling uj ^ q^, this crossover 
is also reflected in the behavior on the plateau in q, Fig. 7: 
for /i > — 0.14ey, the plateau has a negative curvature, 
which can almost be scaled between different dopings, 
while for /x < — 0.14ey, the plateau starts to fill in, ul- 
timately developing a peak at Q. (See also Fig. 3a in 
Ref. 12.) Note that the plateau width collapses in fre- 
quency at both the H- and C-points, while the collapse 
in wave number {qc 0) is only present near the C-point. 
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FIG. 41. Origins of critical cutoffs. Thick solid line = FS; 
thin solid line — Q-shifted FS; dashed lines — Eq. D9, for 
several values of uj. Chemical potential fi = (a) 0, (b) -0.1, 
(c) -0.14, (d) -0.2eV. horizontal arrows indicate uq, vertical 
arrows uj~ . 

From Fig. 41, the critical frequencies (denoted by ar- 
rows) are points where the uj dependence of this length 
changes abruptly, leading to a sharp change in Imx- 
Thus, near the H-point, the plateau width is lo~ (in- 



The parameter B is small, and generally neglected. 
However, it enters into the evaluation of u, so will be 
discussed briefly. The expression for B may be written 
exactly as the uj limit of 



B 



1 



i^k+Q - ^k) ii^k+Q - ^k) 



- ^kf - ' 



(Dll) 



It can be shown that B has a logarithmic correction due 
to the hot spots. The integral can be approximately eval- 
uated by (a) using symmetry to reduce the integral to one 
over an octant of the Brillouin zone containing one hot 
spot, (b) splitting the domain of integration into (i) a 
circle of radius kc about the hot spot, and (ii) the re- 
mainder of the domain, and (c) numerically evaluating 
the integral over domain (ii) while providing an analytic 
approximation to that over (i) . Then the k integral over 
the hot spot circle can be written approximately as 



(1 - ?,(iek)dk 



~2 ["T ~ 3/30 log 



agkc 



(D12) 
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At T = 0, the integral / must then be integrated in 6 
over the wedge where the difference in Fermi functions 
docs not vanish. The integral from outside the hot spot 
circle will eliminate the fcc-dependence, but should not 
affect the log{ui) term. 

It is difficult to directly evaluate the two-dimensional 
principal value integral for B. Instead, it is much sim- 
pler to evaluate Re{x) via Kramers-Kronig transforma- 
tion of Im,(x) and find B by numerical differentiation. 
When this is done, it is found that (a) B is numerically 
very small due to the plateau in Re{x), Fig. 39d, and 
(b) the logarithmic correction is too small to determine 
accurately. 



4. u 

The quartic effective action is 

5 = ^ ^ U2{q,ioJn)(t>{q,ioJn)H-q,-ioJn) 

X(p{qs,iu}3)(t>{q4,i^4), (D13) 

where the prime in the second sum means summing over 
all qi, uJi, such that Yli^i qi = 0, Yli^i = 0, 



■ i(jJ4 {iLVs — Ae){iu!3 + iu!2 — Ae) 



ic02 



{iiOi + Ae)(ia;i -I- iu!2 + Ae) 



1 



{iu>3 — Ae){iijOi + iuj2 + Ae) 
1 

[iuJi + Ae)(zci;3 + iuj2 — Ae) 



where Ae 

{i = 1,2), 
to 



Letting (jji^± = 
and noting that Ui+ = —0J2+ 



{Ui ± Wi+2)/2 

this simplifies 



2U 



where 



W- = (iwi- + Ae) 



)' 



(D18) 
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Thus in Matsubara frequency space, u is largest for 
wi+ = LV2- = 0, so it should indeed be reasonable to 
estimate it in that limit: 



u(zwi,0,0) = U^' 
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d{iu)i)' 



E 



/(gfc+q) - fin) 
iwi + Ae 



(D20) 



ll2{q,i(^n) = -^[1 - Uxo{q,i'^n)], 



(D14) 



n4(9j,«Wn) = E GQ(k,ien)Go{k + qi,ien + iijOl) x 

xGo{k + qi + q2,ien + itoi + iaj2)Go{k - qi,ien - ioJi), 

with u = Il4/Nol3U^. 

Since there is some controversy^^'^^ concerning u, it 
shall be evaluated in detail. Millis''^ showcxl that for free 
electrons (parabolic bands) this expression is in general 
well defined, but diverges when Q is a 'spanning' vector 
of the Fermi surface - in the present case, this would 
correspond to the H- and C-points. Abanov, et al.'^^ 
found a more severe divergence: u diverges for all /x in 
the hot spot regime. The problem lies in the limit of 
external frequencies 0, momenta ^ or Q. Taking 
this limit on the momenta, the expression for u can be 
written as 



In turn, it should be possible to approximate u, Eq. D20, 

by its wi — !■ limit, if this is nonsingular. From Eq. 4, 
Uxo{Q,^) = Blo"^ + iCui + 1 — So- Thus, the analytic 
continuation iuii cu + iS yields 



u(0, 0, 0) = lim 



2BU. 



(D21) 
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The sum over Matsubara frequencies yields 



(D16) 



(D15) 

Due to the plateau in x((5, w), B (Table I) and hence 
u are extremely small. The smallness of u is true only 
in the limit that all external frequencies are small, which 
means that a more complicated expression should be used 
to evaluate u. Moreover, there is an additional problem: 
as found above, B has a correction in In (w), which would 
formally be divergent. Hence, the model is not fully self- 
consistent, and u will be treated as an empirical param- 
eter. The weak logarithmic divergence will be neglected, 
and u approximated by a constant. 



APPENDIX E: INTERLAYER COUPLING 



1. Dispersion of tzi Direct and Staggered Stacking 

Andersen, et al.^^^ demonstrated that the anomalous 
form of interlayer hopping in the cuprates, tz = tzo(cx — 
Cy)2, could be understood by couphng the Cudj-d^ and 
Oj, orbitals to the Cu4s orbitals, which have significant 



32 



interlayer coupling. Here, I provide a simplified calcu- 
lation including only these orbitals, and show how the 
dispersion is modified by staggered stacking of the Cu02 
layers. For uniform stacking (Cu above Cu), the hopping 
matrix becomes 



H 



A 














0^ 





0^ 






-If 



'2t 



(El) 



with Si — smkia/2. Here the first (last) row is for the 
CU(i2_(j2 (Cu4s) orbital, and the middle rows are for the 
Opx and Opy orbitals, with Egz = —'itsz coskzC. In the 
limit As + Efiz » A >> t,tps, the antibonding band 



has dispersion 



E 



A - ^(c. 



2) 



A2(A, + S, 



-{Cx 



(E2) 



so if tsz << As, the interlayer hopping has the form 
tzocoskzcicx - Cy)^, with t^o = -letH^Jsz/^^^l- 
While this form had been suggested earlier^^^ and found 
experimentally for the bilayer splitting in BSCCO^^*, it 
should be noted that it is only approximate, and that, 
at least in YBCO, there is considerable splitting of the 
bilayer bands along the zone diagonal^^^. Nevertheless, 
this form is adequate for the present purposes. 

When successive layers are staggered, the only modi- 
fication to the hopping matrix is in the form of Es{kz), 
which now acquires an in-plane dispersion, 



Es{kz) = -4:tszCoskzc[cos{kx 
which leads to Eq. 65. 



h ky)a/2 + cos {kx — k. 
: cos fczccos kxa/2 cos kya/2, 



2. Estimation of tz from Resistivity Anisotropy 

The dc conductivity can be estimated 



2e2 



^t;f(5(e^-/i)T^, 



(E4) 



i = x,y, z, with Q the unit cell volume, Vi = h^^dej:/dki, 
and the scattering rate. Recent ARPES data suggest 
that, when bilayer splitting is resolved, tj: is relatively 
isotropic over the Fermi surface^^^. Taking indepen- 
dent of fc, the conductivities are given by integrals over 
the Fermi surface. Figure 42a shows a normalized con- 
ductivity ratio. 



ct 



(E5) 



while Fig. 42b shows the resulting normalized interlayer 
hopping izo = t^o \fcja^ which would be required to pro- 
duce a resistivity anisotropy Pzz/ Pxx = 1000. For sim- 
plicity, it is assumed that tzo is small, and azz/'^xx is 
evaluated in the limit tzo — > 0. It can be seen that (a) 
the staggered stacking reduces the conductivity by ap- 
proximately a factor of 20, independent of doping (except 
near the VHS), so (b) assuming the resistivity anisotropy 
is 1000 for optimally doped LSCO, it is estimated that 
tzo/t = 0.11 for staggered stacking; by contrast, if the 
stacking had been uniform, a value oitzo/t — 0.025 would 
have been required. This calculation has recently been 
extended^'^" to the bilayer compound Bi2Sr2CaCu208 
(Bi2212), which has both uniform and staggered stack- 
ing. By comparison with a band structure calculation, 
the approximate values tz/t — 0.38 for the intracell hop- 
ping (corresponding to uniform stacking), and 0.14 for 
the intercell (staggered) hopping were found. 
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FIG. 42. (a) Normalized conductivity ra,tio,azz /(^xx vs 
doping Ef, for uniform (solid line) and staggered stacking 
(long dashed line and short dashed line, (x20)); and (b) re- 
sulting normalized interlayer hopping tzo for staggered (sohd 
line) and uniform stacking (long dashed line and short dashed 
line, (x4.5). 



3. z-Component of Ordering Vector 

Given a finite interlayer hopping tz, the first issue is to 
identify the three-dimensional ordering vector: what Qz 
minimizes the free energy? At mean field level, the initial 
magnetic instability will be associated with the state for 
which the RPA denominator first diverges, i.e., the state 
with the largest value of Rexo{Q,Qz)- (Note that these 
calculations implicitly assume that the two-dimensional 
ground state involves commensurate order at Q.) For 
uniform stacking, a complicated dependence on doping, 
temperature, and tz is found. Figures 43,44 plot xo vs 
chemical potential for T = lOOK, lOK, respectively. The 
shift of the susceptibility peak with doping can readily 
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be understood by comparison with Fig. 1. Both temper- 
ature and interlayer coupUng act to smear out the VHS, 
and in both cases cause the susceptibihty peak to shift to 
smaller chemical potential (lower hole doping), Fig. 43d. 
Note that the peak shifts at different rates for different 
Q^-values, showing that the band is developing a consid- 
erable c-axis dispersion. The fastest shift (short dashed 
line in Fig. 43d, corresponding to Qz = 0) can thus be 
considered as representing a crossover from quasi-two- 
dimensional to fully three dimensional dispersion. 
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FIG. 43. Xo(Q, Qz) at T = lOOif vs chemical potential ^, 
for uniform stacking and Qz = n (a), 7r/2 (b), and (c). 
The various curves correspond to tzo/i = 0.01, 0.02, 0.05, 0.1, 
0.2, and 0.5, with the peak in xo shifting to the right with 
increasing tzQ. Inset (d): position of peak, fimax, vs tzo for Qz 
= TT (solid line), n/2 (long dashed line), and (short dashed 
line) . 
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FIG. 44. xo{Q,Qz) vs chemical potential jj., as in Fig. 43, 
but at T = lOK. 

This dispersive shift of the peak in xo leads to a 
doping dependence of the optimal Qz, as illustrated in 
Fig. 45 for t^o = 0.2t. For large hole doping, near the 
tzO ~ VHS, the susceptibility maximum corresponds to 
Qz — tt/c, while near the susceptibility peak, the spin 
modulation becomes incommensurate (intermediate val- 
ues of Qz have the largest susceptibility). There is a 
rapid evolution of the optimal Qz, and beyond the peak 
regime, over essentially the entire electron-doped regime, 
the optimal Qz is 0. This same pattern is repeated for 
smaller tzo , with only the region of the susceptibility peak 
changing. The results are essentially independent of the 
sign of tz ■ 
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FIG. 45. Xo(Q,Qz) vs chemical potential fj., for uniform 
stacking and t^o = 0.2t, and T = lOK (a), or lOOK (b), 
with Qz/tt = 1 (solid line), 0.75 (long dashed line), 0.5 (short 
dashed line), 0.25 (dotted line), (dot-dashed line). 



4. Calculation of 

a. Uniform Stacking 

Given tz and Qz, the parameter Az of Eq. 63 can be 
evaluated: Ux{Q + qzZ,uj 0) UxiQ + QzZ,0) + 
Az{qz — Qz)"^ ■ The dominant ordering vectors, Qz — tt/c 
and Qz ~ 0, can be analyzed in more detail. For the 
former choice, 



4 



■f k+Q 



ftl rll 

■'k •'k+Q \ 
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with fL = - fj:)/kBT, f'J = -fL{l - 2fj:)/kBT, 

Cz — COS kzC, Sz — sinfczC. For the latter case 
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Figure 46 ( 47a) shows how XoiQ, Qz) varies with Qz 
for tzo — O.li (0.02t), for a number of different dop- 
ings. For the entire electron-doped regime, the peak is at 
Qzm = (Fig. 46b, 47d), crossing over to Qz„i — n/c \n 
the hole doped regime. Away from the peak, the sus- 
ceptibiUty varies as Azq^, with qz = Qz — Qzm, and 
in the electron-doped regime the full variation can be 
approximated by a cosine. The amplitude of the co- 
sine falls to zero as the C-point is approached. In the 
quasi-two-dimensional regime this amplitude scales with 
i^o- Figure 47b, c shows plots of the best parabolic fit to 
A'j, — Azjf? for tz^lt = 0.02 (squares) and 0.1 (triangles). 
For tzviji = 0.1, an alternative is shown, found by fit- 
ting the full susceptibility as a cosine in qz (circles) . The 
good agreement between the two techniques shows that 
this is a reasonable approximation in the electron-doped 
regime {—^.2eV < /x < 0). Near the susceptibility peak, 
the variation is nonsinusoidal, and the parabolic fit leads 
to a large value for A!^. 
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FIG. 46. {a)x<A{Q,Qz) vs for i^o = O.U, and 

T = lOif, and a variety of chemical potentials 
= -0.003559 (solid line), -0.08898 (long dashed line), 
-0.1779 (short dashed line), -0.2669 (dotted line), -0.2847 
(dot-dashed line), -0.2954 (long- long-short-short-short dashed 
line), -0.3025 (long-short-short dashed line), -0.3203 
(dash-dot-dot line), -0.3381 (long-short dashed line), and 
-0.3559 meV (long-short-short-short dashed line). (b) 
= Xo{Q,Qz) — Xo{Q,Qz = 0), where the curves have 
the same meaning as in frame (a). 
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FIG. 47. {&)xo{Q,Qz) vs for t,o = 0.02t, and 
T = lOK, and a variety of chemical potentials fi = -0.003559 
(solid line), -0.08898 (long dashed line), -0.1779 (short 
dashed line), -0.2669 (dotted line), -0.3025 (dot-dashed line), 
-0.3381 (long-long-short-short-short dashed line), -0.3417 
(long-dashed-dotted line), -0.3452 (long-short-short dashed 
line), -0.3488 (long-short-short-short dashed line), and -0.3559 
meV (long-dash-dot-dotted line). (b,c) A'^ = A;^/Uc^ vs /i for 
tzo/t = 0.02 (squares, A'^ x 25) and 0.1 (triangles,circles) . (d) 
Qzm VS fj, for t^o/t — 0.02 (squares) and 0.1 (triangles). 



FIG. 48. {si)A'^ — Az/U(? vs chemical potential ^ for Qz ~ 
(solid lines) or tt (dashed lines) , for a variety of values of t^o 
and T = IQQK. In order of increasing amplitude, the values 
are t^o/t = 0.01, 0.02, 0.05, 0.1, 0.2, and 0.5. (b) Scaling 
of A'j°^ with {tzo/tf. Curves are t^o/t = 0.01 (solid line), 
0.02 (long dashed line), 0.05 (short dashed line), 0.1 (dotted 
line), 0.2 (dot-dashed line), and 0.5 (dot-dot-dashed line), (c) 
Comparison of max{Az) for staggered stacking (solid line) 
and uniform stacking (triangles, xl/20) at t^o/t ~ 0.1. 



b. Staggered Stacktng 

The same calculations can be repeated for the of 
Eq. 65, associated with staggered stacking; Fig. 48a 
shows Az calculated from Eqs. E6, E7 at = (solid 
lines) and tt (dashed lines). The frustration induced by 
staggering of the Cu02 layers is reflected in a strong sup- 
pression of the g^-dependence of x, which leaves a small 
residual contribution quadratic in tzo, Fig. 48b. Since tz 
vanishes at (tt, 0), there is no shift of the susceptibility 
peak with doping. Note the symmetry of the Az values 
between and tt. In fact, xiQz) is closely sinusoidal, 
particularly for small tzo, with maxima either at tt or 
0. Thus, near either the H- or C-points, the maximum 
of X corresponds to Qz — tt. For intermediate dopings, 
Qz = is favored. At two distinct chemical potentials, 
the amplitude of the cosine collapses and changes sign. 
At the crossing points, x is independent of Qz, leading 
formally to Tn 0. Note from Fig. 48c that the sup- 
pression of Az is approximately in the same ratio as that 
of the resistivity, found above. 
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5. Calculation of Tjv 

When there is a finite interlayer hopping t^, Eq.l3 be- 
comes 
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y y 

^ln(^)], (E8) 
ttA I'iD 



where 



5 + Az_q'z^ and T3D = 7r^Az/2Ce'^c^. (A 
small correction to 60 is neglected. Treating the Qz de- 
pendence as a cosine rather than a cutoff quadratic leads 
to qualitatively similar results.) Thus a finite Az always 
cuts off the divergence found in Eq.13, leading to a fi- 
nite Tm whenever there is a zero-temperature Neel state 
(e.g., up to a QCP). It should be noted that the above 
calculation implicitly assumed that T > T^r) ~ Az'. for 
T < T^D the logarithm is cut off and the system be- 
haves like an anisotropic three-dimensional magnet. For 
tzo/t < 0.1, the system is generally in the quasi- two- 
dimensional limit, Fig. 49a. Figure 49b compares the 
mean-field Neel transition with the Neel transition found 
assuming uniform stacking and finite interlayer couplings 
tzo/t = 0.1, 0.02, and 2 x 10"^ [the last found by scahng 
the T^D for tzo/t — 0.02 by the ratio of tl^s]. It is seen 
that T/v ^ as ^ 0, albeit exceedingly slowly. 
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FIG. 49. {a)T-iD vs fi for T = lOK and uniform stack- 
ing with tzo ~ O.lt (circles, xl/25) or 0.02 (squares), or 
staggered stacking with f20=0.02 (triangles), (b) Tn vs x, 
comparing mean field transition (solid line) with interlayer 
coupling models (uniform stacking) assuming tzo/t = 0.1 
(long dashed line), 0.02 (short dashed line), and 2x10"® 
(dot-dashed line) , and the staggered stacking model assuming 
tzo/t ~ 0.1 (dot-dot-dash line). 



FIG. 50. (a) Mean field Tjv vs x assuming paramagnetic 
Ueff (Appendix B). (b) Corresponding Tm vs x, calculated 
using Eq. E8. Squares — staggered stacking with tzo/t = 0.1; 
triangles — uniform stacking with tzo/t — 0.02; solid line and 
circles = data, as in Fig. 31. 



The above calculations are for uniform stacking. For 
staggered stacking is reduced, in approximately the 
same ratio as the resistivities. Hence, the staggered 
stacking with tzo/t = 0.1 should be comparable to uni- 
form stacking with tzo/t = 0.02, as observed. Fig. 49. 
While Tn technically goes to zero for staggered stacking 
near x = —0.0838, the decrease is logarithmic, and in 
practice no more than a weak dip is expected to be ob- 
served (the point with T/v = OK is omitted from the plot 
in Fig. 31). Hence, if tzo is estimated from the resistivity, 
it will be nearly impossible to distinguish uniform from 
staggered stacking via measurements of Tjv. 

In the above calculations, a constant value of A was 
assumed for each doping, as given in Fig. 9. In fact, for 
the electron-doped cuprates, A ~ l/T^-^ for T > T4, 
Fig. 8. This would cause an enhancement of the loga- 
rithmic correction, ^ T^'^, tending to pin Tjv close to 
T^. For the present parameter values, this could reduce 
Tjv by roughly a factor of two, still larger than the ex- 
perimental values. 

A more likely source of the discrepancy is the possible 
temperature dependence of Ueff, Appendix B. The large 
Ueff at half filling arises from lack of screening, in the 
presence of a Mott gap - and is appropriate in analyzing 
the low-T Fermi surfaces found in ARPES. For calculat- 
ing the onset of the Mott gap, the mean field T^r, it is 
more appropriate to use the paramagnetic susceptibility, 
as in Fig 37a. When this is done, considerably smaller 
transition temperatures are found, both at the mean field 
level. Fig. 50a, and when fluctuations and interlayer hop- 
ping are included. Fig. 50b. While the latter are closer 
to the experimental values, no attempt has been made to 
correct Ueff for the short range gap. 
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